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ABSTRACT 

The radiational model used in this study computes the 
planetary albedo and absorption of solar insolation by the 
ocean’s surface and by atmospheric layers for the primitive 
equation model of FNWC. Large-scale cloud parameterization 
in several layers was utilized in these computations. Solar 
insolation disposition was computed from the water-mass and 
the cloud amounts over each gridpoint, Long-wave cooling 
effects were computed using emissivity formulas after 
Sasamori at the earth’s surface and over the same layers, 
and were also found to be dependent upon the cloud para- 
meterization. 

Two forms of the cloud parameterization were tested us- 
ing FNWC data over the oceanic gridpoints for 16 April 1974, 
The objective was to determine the parameterization which 
better verified the radiational balance as a function of 
latitude when compared with satellite climatology at the 
top of the atmosphere for the same general data-period in 
1969 (after Raschke, et al, 1973), The better verification 
resulted with the smaller cloud parameterization values. 
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I 



INTRODUCTION 



This thesis is a study of a radiative heating parameteri- 
zation for use in the Fleet Numerical Weather Central (FNWC) 
prediction system. The study has as a primary objective, the 
analysis of the radiational and heat balance of the ocean- 
atmosphere system utilizing FNWC gridded data-fields at con- 
stant pressure levels for 16 April 1974. The gridpoint data 
were selected along four oceanic meridians, three in the 
Pacific and one in the Atlantic, with most of the gridpoints 
located in the Northern Hemisphere (Fig. 1). 

The specification of amounts of clouds in two designated 
layers has the most influence on the radiative-model disposi- 
tions (short- and long-wave). The initial specification of 
the fractional amounts of CL(1) and CL(2) are based on large- 
scale formulations developed by Smagorinsky (1960) and used 
in a similar study based on the data day 16 January 1974 by 
Spaeth (1975). 

The governing equations of the radiational transfer of 
the model have been derived by Martin (1972, 1974), who modi- 

fied both the solar and terrestrial radiation transfers to 
respond to the presence of clouds. These clouds were of 
specified amounts in the two layers, one of which is a mid- 
level cloud and the other a low-level cloud. The radiational 
model in use here has similarities to those in use in UCLA 
and NCAR General Circulation models. 
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The physical description of the radiative model may be 
applied to any scale of analysis for which there is adequate 
resolution of the temperature and moisture data in the verti- 
cal. In the horizontal, the reliability of the data used 
here is consistent with that of the FNWC analysis to grid- 
points, and is typically reported to the nearest tenth of a 
degree with regard to temperature and dew point. The radia- 
tional computations made here are applied to FNWC gridpoints 
and are designed to make a one-hour forward-time step appli- 
cable to the FNWC primitive equation forecast model, with 
special adaptations to their a-levels. 

The Smagorinsky (1960) cl oud- sp ec if i ca t ion of CL(1) and 
CL(2) when used with the model of this study resulted in 
global albedo estimates that were too high as compared to 
the satellite data of Raschke (1973), thus leading to under- 
estimates of the net radiative balance at the top of the 
atmosphere and at the ocean surface. This problem suggested 
the reduced cloud-amounts model formulated by Spaeth (1975) 
and used in this study. 

The comparative results afforded by the 2/3-CL parame- 
terization gave reasonably close agreement with the radia- 
tive climatology of Raschke et al (1973), for the NIMBUS III 
period 1-15 May 1969. Oceanic and tropospheric balances 
were computed for both sets of cloud conditions (full-CL 
and 2/3-CL) for 16 April 1974. It was also necessary to 
augment the purely radiative model by including a turbulent 
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boundary-layer model for sensible and latent heat transports 
over the ocean as adapted from Kaitala (1974). 

The. radiation package results computed with the 2/3-CL 
cloud model compared very favorably with satellite clima- 
tology, though only one data day (16 April 1974) was used as 
representative of mid-April soundings over the ocean, both 
for the model and the Raschke climatology. However, the 
latter was based on the data period 1-15 May 1969. 

The major effect noted by the use of the 2/3-CL as con- 
trasted with full-CL parameterization was the reduction in 
global albedo and the resulting greater contribution to the 
surface net heating rates. This result is in general agree- 
ment with recent observations by Van der Haar and Hanson 
(1969), Van der Haar and Oort (1973) and Raschke et al (1973). 
However, the major climatological discrepancies resulting 
from the specification of layered cloud amounts in the 
tropics in a manner identical to that used poleward of the 
subtropics still seems to give somewhat excessive albedo. 

It is felt that the geometry of solar radiation streams im- 
pinging upon subgrid sized cumulus cells in the tropics 
should be remodeled so as to divert a greater percentage of 
solar radiation downward, as compared with the characteris- 
tic cloud-layering and resultant reflective effects which 
seems permissible in middle and high latitudes. 
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II 



DATA PREPARATION 



A. INITIAL DATA FIELDS 

The temperature and humidity data used in this study were 
arranged in the form of soundings taken along four oceanic 
meridians (Fig* 1) of the Fleet Numerical Weather Central 
(FNWC) Northern Hemisphere analyses on 16 April 1974* Compu- 
tation of radiational dispositions and of other heat budget 
terms are made at these gridpoints* Oceanic locations for 
these computations were chosen because: 

p 

1. The constant O surfaces (where a - — ) of the FNWC 

TT 

primitive equation system are close to being constant pres- 
sure levels • 

2. The maritime-area heating rate computations are like- 
ly to be representative of the month of April 1974 as con- 
trasted with computations made for 16 April 1974 over a 
corresponding set of gridpoints over land. 

The three meridians (and their respective number of sound- 
ings) selected over the Pacific Ocean were located at 125W 
(25 soundings), 170W (25 soundings) and 145E (17 soundings). 
The Atlantic Ocean meridian was 35W (26 soundings). This 
method of selecting "soundings" along the indicated meridians 
of the FNWC polar stereographic map, made it unnecessary to 
employ spatial interpolation between original data gridpoints 
along the meridians. Data along line 3 in the Pacific was 
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Figure 1. FNWC polar stereographic grid and meridians (lines 
1, 2, 3, and 4) selected for study. The longitudes A are 
shown for each meridian as well as the extent considered of 
each meridian. 
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not extended southward of gridpoint (9,55) because they fell 
over land masses (New Guinea and Northern Australia) where 
the surface temperatures and other sounding features were 
unrepresentative of the oceanic values* 

The gridpoint soundings in the form of Table 1(a) were 
taken from the original FNWC 63-by-63 surface analysis, as 
well as of the vertical distribution of T(P) and of dew 
point depression at five other standard pressure levels up 
to and including 400 millibars. These original soundings 
were then modeled into a radiative sounding (Table 1(b)) cor- 
responding to FNWC primitive equation prediction levels up 
to 100 mb (see Fig. 2). 

Data from 0000GMT, 16 April 1974 were used for the Paci- 
fic soundings while 1200GMT, 16 April 1974 data were used 

for the Atlantic soundings. These times were used so that 
the analysis times would correspond as closely as possible 
to local solar noon in each of the indicated areas. 

At each gridpoint selected, the original humidity data 
was given in the form of five dew point depressions for the 
analysis levels from 925mb to 400rab. Since the standard in- 
strument level vapor pressure ( e a ^ r ) was missing from the 

original FNWC analysis at all gridpoints, e . was approxi- 

a x r 

mated using the FNWC field of e^, a computed value of vapor 
pressure at about 20 meters above mean sea level (MSL) in 
the turbulent boundary layer. FNWC values of and e^ over 

the ocean had been made available using an operational 
planetary boundary layer model detailed by Kaitala (1974). 



23 



Table 1(a) , Example of typical FNWC sounding for gridpoint 
(1^1). The humidity parameters between 925,.., 
400mb are dew point depressions, but at the surface and top 
of the constant flux layer (9999*) are vapor pressures in mb. 



Pressure 


(mb) 




t (°e) 


Humidity Parameters 


Surface , 


1000 mb 




25.6 


e = 27.4 mb 

X 


925 






18.3 


1.2 C 


850 






14.6 


1.0 C 


700 






5.6 


3.4 C 


500 






-11.0 


2.8 C 


400 






-22.7 


• 1.4 C 


300 






-37.6 




250 






-46.5 




200 






-56.3 




150 






-67.3 




100 






-80.6 




9999* 




Tx 


= 23.5 




* 

Code 9999 


indicates 


data 


taken from 


the top of the constant 


flux level 


(CFL) of 


the 


FNWC initial 


data program. 



Table 1(b) , Example of the corresponding radiative sounding 

with mixing ratio listed at odd k-levels (Fig. 2). 
Additionally, water vapor and CO^ absorber masses, cloud 
amounts, CL(1) and CL(2) (after Smagorinsky (I960)) as these 
parameters are modeled in the radiative theory. 



r essure 
(mb) 


T(°C) 


Mixing 

Ratio 

g/kg 


Absorber Masses Smagorinsky 

water vapor CO^ Cloud Amounts 

(gm/cm^) ( c ni/cm^) CL(1) CL(2) 


1000 


25.60 


17.10 








900 


16.67 


12.56 






l.ooo, cl 2 


800 


11.84 




2.35 


45.53 




700 


5.60 


6.43 








600 


-1.57 




3.36 


83.53 




500 


-11.00 


2.68 






0.900, CI^ 


400 


-22.70 


1.41 


3.69 


113.35 




300 


-37.60 


0.68 








200 


-56.30 




3.75 


133.99 




100 


-80.60 


0.04 


3.75 


138 . 67 




0 


-80.60 


.00 


3.75 


143.35 
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Level - 
Ideniifers 

0 



U(0), C(0) 



O-level 



T 



Layer (0,2) 



0 = 0 



1 



T.q 



0.1 



U(2), C(2 ) 



0 = 0.2 




3- 



T.q 



0.3 





U(10), C(10) 



0 = 1.0 



Figure 2. Five-layer radiative sounding used in this study. 
Levels are identified by their values on the k-scale, while 
layers are identified by their level boundary indicies in 
parentheses, e.g. (8,10). Pressure-scaled water vapor and CG 2 
mass increments dU and AC , respectively are integrated rel- 
ative to the surface and the resulting U and C are carried in 
the model at even levels. The temperature T is retained at all 
levels. Amounts of clouds CL(.l) and CL(2) in the layers shown 
have been parameterized for consideration of their radiative 
effects . 
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To perform radiative calculations, it is necessary to 
have water vapor and CC^ absorber masses and cloud amounts, 
CL(1) and CL(2) at certain required k-level boundaries 
(Fig. 2). All soundings in this study start at sea level 
with the approximation of surface pressure tt = lOOOmb. 
Therefore, the k-levels correspond approximately to the 
FNWC levels P fc = 1000., 900., 800., ...., 200., 100., O.Omb 

corresponding respectively to = 1.0, 0.9, 0.1, 0.0. 

B. INTERPOLATIVE PROCESSING TO k-LEVELS IN RADIATIVE 
SOUNDINGS 

1 • Temperature Profiles 

The gridpoint temperatures were listed at each manda- 
tory level of Table 1(a) between 1000 o , ..., lOOmb. The tem- 

perature was assumed to be isothermal from lOOmb to p=0.0mb. 
The temperature T^q was set equal to the FNWC listed sea-sur- 
face temperature. The radiative sounding temperatures for 
the remaining k-levels were obtained from either their corres- 
ponding listed temperature-level or by a three-point Lagran- 
gian interpolation scheme (Eq. 2-1, Spaeth, 1975), to level 
k when the listed FNWC temperature-profile did not include 

the value T, . 

k 

2 . Moisture Profiles 

The moisture from the original FNWC soundings (Table 

1(a)) were converted into mixing ratios at each of the 

original sounding levels. The near-surface vapor pressure 

(taken as e ) was used to calculate the mixing ratio at k=10 
x 
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as described by Spaeth (1975): 



q 10 “ 621.97 (e x /1000) (2-1) 

The remaining ratios q were calculated from the original 
FNWC sounding (Table 1(a)) dew point depression data after 
Spaeth (1975, Eq. 2-3). The computed q-values were subse- 
quently interpolated to k-levels using the previously men- 
tioned three-point Lagrangian procedure applied to successive 
q-values in the original sounding. Resulting values of the 
interpolated q-values are shown for the case of the radia- 
tive sounding at gridpoint (1,1) (Table 1(b)). 

Due to the fact that most radiosonde humidity data 
for p 300mb is either unreliable or not available, a power- 
law extrapolation formula (Spaeth, 1975) 

is /k s A 

q 5 = 500^ (2-2) 

was used to obtain q-values at k=*3,2,l. 

This extrapolative procedure for obtaining q-values 
at high levels was first suggested by Smith (1966), utilizing 
the least squares estimation of X according to 

6 

2 y 

X - 

E x 
i 

<l ± 

with y, = log and x. * log(P./P ). The Smith method was 

l q ^ l id 

then further tested by Spaeth (1975) for the vertical scale 



. x . 

l l 



2 R yx 






£ x . 

i 



(2-3) 
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of the FNWC initial data soundings. The resulting correla- 
tion coefficient R was generally found to lie in the range 
.95 to .99. Such high correlations indicate that the pro- 
cedure described by Spaeth (1975) for determing upper atmos- 
phere q-values has realistic estimation value for the 
vertical scale involved in the radiative sounding. 

3 • Pressure-Scaled Absorber Masses 

The pressure-scaled water vapor absorber mass Au in a 
layer (Fig. 2) was calculated for the five odd numbered k- 
levels by using the computed mixing ratio values (Eq. 2-7, 
Spaeth, (1975)). The equation for computing the integrated 
pressure-scaled water vapor mass is given by the algorithm 
described by Eq . 2-8 of Spaeth. 

A similar algorithm was used in the computation of 
the carbon-dioxide scaled absorber mass (Spaeth, 1975, Eq. 
2-10,11). To be noted here is the fact that CO^ absorber 
masses have been stated in terms of N.T.P. "pressure-scaled" 
cms, since the CO^ absorption coefficients are generally 
stated in terms of this CO^ mass unit. 

C. CLOUD PARAMETERIZATION 

The relative humidities, and thus the saturation vapor 
pressure at levels k=5 and k=9, are used in the calculations 
of the fractional cloud cover CL^ and CL^ in layers (4,6) 
and (8,9) respectively, as depicted in Fig. 2. The fractional 
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cloud amounts for these two layers were parameterized using 
the following equations (after Smagorinsky, 1960): 

CL ( 1 ) = 2.0 (RH ( 5 ) ) - 0.7 (2-4(a)) 

CL ( 2 ) « 3.33 (RH (9 ) ) - 2.0 (2-4(b)) 

Smagorinsky 1 s parameterization of CL(1) and CL(2) permits 
cloudliness fractions of 1.30 and 1.333 respectively with 
RH=1.0. Values of CL greater than 1.0 were considered 
initially by Smagorinsky (1960) to suggest the amount of 
supersaturation which accompanies precipitation. 

All operational radiative models in use at the present 
time limit CL 1.0. In the radiative study developed here 
no opportunity is afforded to deduce precipitation-rates 
and/or super saturation amounts. Hence, the Smagorinsky for- 
mulations of CL^ and CL 2 were reduced in each case by the 
multiplicative factor 2/3 so that neither fractional cloud 
cover could exceed unity. Thus alternative cloud-cover 
fractions 



CL ’l ■ 


2/3 


(CL 1 ) 


(2-5 (a) ) 


CL, 2 " 


2/3 


(cl 2 ) 


(2-5 (b) ) 


cons id er ed 


for 


purposes of 


estimating the shift in the 



global radiative balance as a result of altering the cloud 
computations in the radiation model. 

The 1/3 reduction inherent in CL f , CL 1 ^ relative to 
Smagorinsky 1 s CL^, CL^ were initiated in this study in an 
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attempt to tune cloud amounts to give radiation results in 
closer agreement with recent satellite climatology of 
Raschke et al, (1973). This satellite climatology suggests 
the use of smaller cloud amounts than that specified by 
Smagorinsky . Fractional cloud amounts were considered to 
be functions of large scale effects only. Therefore, small 
scale convective activity, seasonal conditions (except for 
the soundings) and latitudinal effects were not considered 
in specifying the reduction factor of Eq . 2-5. The reduced 

cloud-parameterization was introduced here for estimating 
large scale radiational effects only. 

D. CLOUD-AREA COVERAGES 

Since Eqs. 2-4a,b or 2-5a,b gave the fractional-coverage 
of the gridpoint area by the appropriate cloud-type, the 
gridpoint area may be thought of as broken into random frac- 
tional segments of size 

W(0,0) = (1“CL 1 ) (1-CL 2 ) ( 2-6a ) 

wherein there is a combination of clear-over-clear segments 
in the layers. Similarly, the gridpoint area has the frac- 
tional area of cloud coverage 

W(l,l) = CL * CL 2 ( 2-6b) 

of upper-cloud amount CL 2 overlying lower-cloud amount CL 2 . 
Likewise the combinations of cloudy over clear and clear 
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over cloudy , by layers, may be visualized as occurring with 
the respective weights 

W(1,0) = CL 1 * (1-CL 2 ) ( 2-6c ) 

and 

W(0,1) = (1-CL 1 ) * CL 2 ( 2- 6d ) 

Regardless of whether the full-cloud amounts CL^, CL 2 of 
Smagorinsky (Eq. 2-4), or the 2/3-CL amounts of Eq. 2-5 were 
utilized, it was useful to carry, for radiation computations 
of each sounding, the relative weights or fractions of the 
gridpoint area exposed to the specified cloud-layer combina- 
tions. Henceforth, the symbols denoted by W(0,0), W(l,l), 
W(1,0), W(0,1) and as given by Eqs c 2- 6 ( a , b , c , d ) , suggest 
overcast (1) or clear (0) cloud amounts in the indicated 
layers (Fig. 2), the first index 1 or 0 referring to layer 
CL^, and the second to CL 2 * 

The usefulness of this computational device will be 
clarified in Sections III and IV, where the procedures for 
the terrestrial and solar radiational computations are dis- 
cussed and the results are summarized over the set of 
soundings . 
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Ill . TERRESTRIAL RADIATION 



A. THEORETICAL AND EMPIRICAL BASIS 

Empirical formulas were developed by Sasamori (1968) for 
flux emissivities in the atmosphere associated with computa- 
tions for the heat balance requirements of the NCAR General 
Circulation Model. Sasamori derived the empirical emissivity 
formulas for water vapor and CO^ by comparison with the 
theoretical values built into the Yamamoto Radiation Chart 
(1952). This chart has proved to be quite accurate for com- 
putational checking of the Sasamori emissivities and was used 
by Spaeth (1975) and Warner (1974) as a systematic guide for 
integration of the radiative transfer formulas developed by 
Martin (1972, 1974), who adapted the Sasamori formulas to 

the particular layers of interest in the gridpoint computa- 
tions of the FNWC primitive equation model (Fig. 2). 

The essential long-wave flux formulas required for use in 
the FNWC heating package are the following: 

net IR flux at earth, k=10 
net IR flux at level k^ 
net IR flux at level k=2 

net IR flux divergence in the layer (6,10) 
net IR flux divergence in the layer (2,6) 



F 10* " 
V = 

F * =. 

2 

F610 = 
F26 = 



To get the IR flux divergences in the layers (6,10) and 
(2,6), the differences Fg*-F^Q* an< * F 2* -F 6* must be computed 
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The detailed scheme for making such computations with various 
combinations of cloud cover CL(1) and CL(2) is similar to the 
technique developed by Martin (1974) and as reproduced in 
detail by Spaeth (1975). 

In order to make IR net-flux calculations along the path 
of integration, there must be a physically sound representa- 
tion of the emissivity (e ) as a function of both water 
vapor and CO^ absorber masses in layers along the sounding. 

For a further discussion of the emissivity formulas used in 
the quadrature scheme, refer to Spaeth 1 s Appendix A (1975). 

B. NET FLUX FORMULATION 
1 . At Level k s 10 

The radiative sounding as depicted in Table 1(b), was 
computed as the combination of parameters U(k,10), C(k,10) 

and T for each required level, k=*l 0 , 8 , . • • , 2 , 1 , 0 . Cloud para- 

K 

meters CL(1) and CL(2) are also listed at each gridpoint and 
in general are non-zero. The grid area was then considered 
to be composed of areal fractions (weights) defined in 
Eqs. 2-6(a,...,d) and denoted by the symbols W(0,0), W(l,l), 
W(1,0) , W ( 0 , 1 ) . 

The composite net flux (CL^, CL 2 ) at level k=10 

at each gridpoint is then constructed by using the appropriate 
weight factor to multiply the ref er ence net flux F^* compu- 
tations defined for the special cloud-cover cases 

F 10 *(0,0), F 10 *(1,0), F 10 *(0,1), F 10 *(l,l) 



34 



so that 



F 10* (CL 1’ CL 2 )=W(0 ’ 0)F 10* (0j0)+W(1,0)F 10* (1,0) 

+W(0,1)F 10 *(0,1)+W(1,1)F 10 *(1,1) (3-1) 

Spaeth (1975) has listed these reference net flux formula- 
tions in his Eqs. 3-6, 7,8. Using the definitions of W(0,0), 

W(1,0), W(0,1), W ( 1 , 1 ) , (CL^ , CL^) can be shown to assume 

the form 

F 10* (CL 1’ CL 2 ) = U-CL(2) ] { (B 1Q -B 6 )-.5[e wc (8,10) (B 10 -B g ) 

+ (e wc (8,10)+ e wc (6,10))(B 8 -B 6 )]} 

+ (1-CL(2)) (l-CL(l)){B 6 -.5[(e wc (6,10) 

+ e (4,10)) (B -B )+(e (4,10) (3-2) 

w c oh w c 

+ e wc (2,10)) (B 4 -B 2 )+(e wc (2,l0) 

+ e wc (l,0)) (B 2 -B 1 )+ e wc ((0,10),T 1 )*B 1 ]} 

+ CL(2){(B 10 -B 9 ) [l-.5e wc (9,10) ]} . 

Here 

B k = 1.170403 x 10 -7 T k 4 (3-3) 

is the Stef an-Boltzmann blackbody flux in langlies per day. 

Further, e (U. , C. ,10) is the combined water-vapor 
wc k k 

and CO^ emissivity along the path from level 10 to level k. 
This emissivity is considered by Sasamori to be temperature 
independent for T >L 210K, whereas represents the 
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temperature dependent emissivity applicable for T < 210K 
(see pp. 136-137, Spaeth, 1975). 

The reference net fluxes °f Eq . 3-1 are asso- 

ciated with (1) clear skies in both layers, (2) overcast in 
the upper layer only, (3) overcast in the lower layer only 
and (4) overcast in both layers respectively. 

2 . Net Flux F 6 _* 

The formula for F *(CL-,CL 0 ) has been developed by 

o 12 

Spaeth (Eq. 3-10, 1975) in a manner analogous to the deriva- 

tion of the weighted F^q*. The result is reproduced after 
Spaeth (1975, see pp. 43). 

F 6 * = tl-CL(l)]{B 8 -.5te wc (6,8)(B 8 -B 6 )+e wc (4,6)(B 6 -B 4 ) 

+ (e wc (4,6)+ £wc (2,6))(B 4 -B 2 ) + (£ wc (2,6) 

+e wc (l,6)) (B 2 -B 1 )+ e wc ((0,6) ,T 1 )*B 1 ] } 

+ ( 1-CL ( 1) ) ( 1-CL (2) ) { (B 10 -Bg) [l.-.5(e wc (6,8) 
+£ wc (6,10))]}+CL(1){(B 8 -B 6 )* (3-4) 

[l-.5e wc (6,8) ] }+CL ( 1 ) (1-CL(2)){ (B 10 -Bg)* 
[l-.5e wc (6,8)+e wc (6,10))]} . 

3. Net Flux F * 

The net flux at level k=2 can be calculated in a 

similar fashion to F_ * and F *, with the formulation of the 

1 U o 

individual reference net fluxes as previously described. 
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The formula for F 2 *(CL^,CL 2 ) is then developed analogous to 

that for F * (CL- , CL 0 ) and F * (CL , CL 0 ) , 

1U ± Z O 1 z 

F 2*~ tl-CL(l) ]{Bg-.5[e wc (2,4) (B^B 2 )+ (e^ ( 2 , 4 ) 

+ e wc (2,6))(B 6 -B 4 )+( £wc (2,6)+e wc (2,8))* 

(B g -B 6 )+e wc (l, 2) (B 2 -B 1 )+e wc ((0,2) ,T 1 )*B 1 ] } (3-5) 

+(1-CL(1))(1-CL(2))*{(B 10 -B 8 ) [l-.5(e wc (2,8) 

+ e wc (2,10))]}+CL(l){B 4 -.5[e wc (2,A) (B 4 -B 2 ) 

+ e wc (l,2) (B 2 -B 1 ) + e wc ((0,2) ,T 1 )*B 1 ]} . 

A new parameter, the total outgoing long-wave radia- 
tion to space (FF2) can be readily defined from the expres- 
sion for F 2 * by setting to zero the terms of Eq. 3-5 
representing the downward flux through level k=2. These 
terms are just those denoted by the symbols & wc ( (0 , 2) , T^) *B^ 
and (1 , 2 ) (B 2 “B^) . If the long-wave radiative model, in- 

cluding its cloud parameterization, is realistic, the modeled 
outgoing flux to space (FF2) should compare reasonably with 
that observed by NIMBUS III satellite during the same period, 
after Raschke et al (1973). The most nearly comparable 
satellite climatological period was 1-15 May 1969. These 
satellite observations were obtained from the NIMBUS III 
satellite atlas of Raschke et al (1973) for the May dates 
indicated and used for comparison with the model values of 
FF2, point for point (at 5°latitude intervals) along the 
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same four meridians of Fig. 1. Also the la t itudinally- 
distributed mean satellite fluxes were compared with the 
corresponding averaged values of FF2. The results showing 
the comparisons of the model computations of FF2 with the 
corresponding satellite observations are listed by latitude 
in Table V. 

C. APPLICATIONS TO HEAT BALANCE COMPUTATIONS 
1 • At Upper and Lower Boundaries 

In order to compute the heat balance at the top of 
the ear th- tropospheric system the composite Y and the 
total insolation absorbed below level k=*2 at each gridpoint 
are required. Moreover a radiative balance is computed at 
the surface for each gridpoint using the earth 1 s absorbed 
insolation (Section IV) and the composite net flux loss F^^* 
In order to consider the possibility of a heat 
balance at the earth’s surface, rather than merely a radia- 
tive balance, it is necessary to add a term representing com 
bined sensible and latent heat transfer across the air-sea 
interface (Section V). This latter transfer rate at each 
gridpoint was adapted from the FNWC primitive equation model 
after the discussion of Kaitala (1974) and based upon the 
FNWC data of 16 April 1974. 

2 • Intermediate Levels 

Consideration of the radiation balance in the atmos- 
phere requires computation of the long-wave cooling effects 
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caused by the flux divergences F26 and F610. F26 and F610 

are defined as 



F 2 6 = F 2 * - F 6 * 


(3-6) 


F610 = F * - F * 
0 1 u 


(3-7) 



where the symbols (2,6) and (6,10) indicate the layer 
boundaries involved in the balance considerations for the 
indicated layer. 

In order to compute F26 and F610 at each gridpoint, 

complete listings of IR net fluxes F^*, F^*, ^ 2 * 9 associated 

with each set of reference cloud amounts (0,0), (1,0), (0,1), 

(1,1) have been computed at each gridpoint together with the 

weighted set of fluxes F_ * ( CL- , CL 0 ) , F , * ( CL - , CL 0 ) , 

1 U I/O 1 Z 

(CL^ , CL^) • These cloud-weighted values have been con- 
structed using the weighting scheme of Eqs. 2-6, 3-1, 3-4 and 

3-5. A sample gridpoint printout of the net fluxes has been 
included in Table II. 

The determination of the heat balance by layers has 
been deferred to Section VI, where the results are displayed 
in cross-sectional form for the 2/3-CL case. 

D. STATISTICAL RESULTS AND COMPARISONS 
1. Net Flux F *(0,0) 

A statistical test of the F^q*( 0,0) numerical results 
was conducted by utilizing a linear regression program for 
relating the predictand F^q*( 0,0) against the simultaneous 
values of and X 2 =B 10 V ^ as P red * ctors > using all 93 
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gridpoint data as samples. The regression program was based 
on the BMD02R in the Biomedical set of programs (Dixon, 1973). 
The result was cast in the form of the well-known Brunt net 
flux equation (1932) : 



where 



F 10* - B 10 (a + b/a) 



10 



St * T 




(3-9) 



is the S t ef an-Bo 1 tzmann blackbody flux at the surface tempera- 
ture and e is the surface vapor pressure in mb. The 

multiple correlation coefficient R ^( F 20* ^ X 1 ’ ^2^ W3S higher 
than 0.99 in the 16 April case, as was also true for the data 
periods of Spaeth (16 January 1974), and Beahan (16 July 1974), 
(Table III) . 

Since F^q* may be defined as 

F * = R - F (3-10) 

10 10 d v ' 

Eq . 3-10 can then be solved for F^ 

F d 3 B 10 ((l-a)-b/e} , (3-11) 



for the clear - sky case. It must be recalled that the compu- 
tation of F, included both the effects of water-vapor and of 
a 

CO^. Since there is no direct relationship of /e upon the 
quasi-constant absorber mass of the radiative effect of 

the latter constituent in Eq . 3-10 must be included in the 

term involving (1-a) B^q. 
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CL (1 ) = 0.0 CL ( 2 ) = 0.0 WT - 0.0 F * = .1635 F * = .2645 F * = .3483 ly/min 
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A cross-seasonal comparison of the coefficients 1-a 
and -b of Eq. 3-11 is presented in Table III, comparing the 
clear-sky results of for 16 April 1974 with those of Spaeth 
and Beahan at intervals of three months earlier and later, 
r espec t i vely . 

The high multiple correlation coefficients indicate 
that an oceanic version of the Brunt downward flux equation 
is valid for each tested season. depends most strongly 

upon B 1q while only slightly on B^^/e. This is an apparent 
result of the nearly constant relative humidity over the 
large oceanic regimes. Moreover, the weak dependence of F^ 
upon B^g/e is most evident at the data-time 16 July 1974, 
when the coefficient -b is smallest. This characteristic is 
due primarily to the relatively weak gradient of /e in the 
radiative soundings aloft. 

16 January 1974 16 April 1974 16 July 1974 



1-a 


.6436 


.6557 


.8129 


-b 


.0224 


. 0209 


- .0094 




.4500 


.4527 


.5471 


r m 


.9955 


. 9969 


. 9950 



Table III. Cross-seasonal comparisons of the coefficients 
1-a and -b from Eq . 3-11. The mean values F^ 

from the Brunt equation and the resulting multiple 
correlation coefficients have been included. 

2. Modification of F ^ q * f ° r Cloudiness C L ^ , C L ^ 

A second statistical test was performed based on a 
possible relationship between the cloudy-sky and the clear-sky 



cases of F 1Q *. The ratio Y = F 1Q *(CL , CL 2 ) / F Q * ( 0 , 0 ) was 
used as the predictand in this case. The predictor CL was 
the total opaque cloud cover at each gridpoint, defined as 
(after Quinn, 1971 and Spaeth, 1975): 

CL = CL ( 1 ) + CL ( 2 ) (l-CL(l)) . (3-12) 

This regression was developed utilizing (CL^ , CL^) in 

the numerator, computed with the amounts CL(1) and CL(2) as 
calculated from Smagor insky 1 s formulation, and then utiliz- 
ing F 1 q * (CL 1 -pCL 1 2 ) computed from the 2/3-CL parameterization. 
The general form of the regression formula required is shown 
in Eq. 3-13 below, with CL given by Eq. 3-12 above. 

F 1q *(CL) = F 10 *(0,0) [ 1- d ( CL ) ] . (3-13) 

A cross-seasonal comparison of the coefficients d and other 
statistics for the 2/3-CL and full-CL cases using the results 
of this study and those of Spaeth and Beahan, is presented 
in Table IV. There seems to be little clear-cut variation in 
the statistics of Table IV between the results of adjacent 
seasons, and little statistical preference between the full- 
CL and the 2/3-CL cases. 

The high correlation coefficients of Table IV for 
all seasons considered, indicate the general capability of 
the IR radiative model employed in this study to account for 
the effective net radiation at the surface. 
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16 January 1974 16 April 197A 16 July 1974 

2/3-CL Full-CL 2/3-CL Full-CL 2/3-CL Full-CL 

d .75 .75 .73 .75 .79 .81 

R m .9936 .9916 .9945 .9942 .9920 .9916 

Table IV. Cross-seasonal comparison of the coefficient d 
from Eq . 3-13 and the resulting multiple corre- 
lation coefficients for both the 2/3-CL and the full-CL 
cases . 



The mean values of 
ings for 16 April 1974, as 
were as follows: 

Full-CL Statistics 
F 1 q * (CL) = .0639 ly min' 1 
F 1() * (0,0) = .1515 ly min" 1 
CL = .7733 



Fio* over the 93 gridpoint sound- 
applicable in Tables III and IV 

2/3-CL Statistics 
F 1Q * (CL) = .0868 ly min" 1 
F 1Q * (0,0) = .1515 ly min" 1 
CL = .5888 



The superior bar symbol denotes the sample mean. 

These results indicate that the surface net flux is decreased 
by the respective ratios .578 with CL = .7733, and by .427 
with CL = .5888. These results indicate that the downward 

flux ( CL ) is 57.8 percent larger at the surface than 

Fj(0,0) when full-CL amounts are used in the cloud parameteri- 
zation model. In the 2/3-CL parameterization the mean down- 
ward flux at the surface F^(CL’) is increased by only 42.7 
percent when compared to F^(0,0). 

E. COMPARISONS OF FF 2 WITH SATELLITE CLIMATOLOGY 



Comparison was made of computed-model values of FF2 with 
satellite measurements of total long-wave flux to space 
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(after Raschke et al, 1973) for the NIMBUS III period of 
1-15 May 1969* FF2 is used here instead of F^* since the 
downward flux at k=2 is not measured by the satellite. The 
computation of FF2 was previously discussed. FF2 for 16 
April 1974 was computed for both full-cloud and 2/3-CL cases 
at all 93 gridpoints. These FF2 values were then inter- 
polated to whole multiples of 5 degrees of latitude along 
each of the four meridians and the results were then averaged 
to get a mean latitudinal distribution value for each 5 
degree latitude interval considered in this study. 

Table V shows both of the la t i t ud inal ly distributed FF2 
cloud-model values compared with those extracted from 
Raschke et al (1973). Raschke’s results were obtained by 
interpolation from charts using the same oceanic meridians 
as those used in this study. La ti tudinally distributed total 
opaque cloud cover CL by Eq. 3-12 are also listed for both 
the full and 2/3-CL cases at each latitude. In the bottom 
line of each column in Table V is listed the cosine weighted 
mean of each set of the column values for the listed latitude 
range (see Eq . 7-2 for the equation defining "weighted 

average") . 

It should be noted that the Raschke results do not corres- 
pond to specifically known CL values. Both cloud cases give 
FF2 results which are reasonably close to the values reported 
by Raschke, especially the means and between 5N-25N. Neither 
the full-CL nor the 2/3-CL model can be conclusively selected 
as preferable, based on Table V results alone. The limitations 
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of the comparisons made here are obvious, when it is re- 
called that between latitudes 20S-5S and between 60N-65N 
there are fewer than four meridional lines available for 
computing the listed zonal values in Table V. For all other 
latitudinal average values, four meridional lines were used 
in the averaging. 

F. COMPARISONS OF CROSS-SEASONAL RESULTS OF F 1Q * 

Table VI depicts the la ti tudinally-dis tr ibuted values of 
F 1q * and of CL at the earth’s surface obtained using the 2/3- 
CL parameterization. A seasonal comparison of the model com- 
puted F^q* is made with data obtained from Spaeth’ (1975) for 
the 16 January 1974 and from Beahan (1975) for 16 July 1974 
cases. It is clearly shown that F^q*(CL ! ) is a decreasing 
function of cloud cover. There is a clear-cut tendency in 
each season for a maximum value of F^q* to be located in the 
subtropics (latitudes 15N-25N). Also there is evidence of a 
high latitude (55N-60N) minimum F^q* associated with a con- 
centration of maximum cloud cover CL 1 . The most outstanding 
variation is the transition in the Southern Hemisphere lati- 
tudes (20S-10S), which has large cloud cover in both the 
January-Ap ril monsoon period but only small CL T in the July 
data period. The same relative conclusions resulted when 
the full-CL parameterization was examined for cross-season 
comparisons, but the obtained results are not included in 
Table VI. 
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Lat. 


(Full-CL) 


FF 2 

(2/3-CL) 


F2 (RAS) 
(Ras chke ) 


C lo ud- 
Full-CL 


Amount s 
2/3-CL 


20S 


' .2581 


. 2852 


. 41 


.985 


. 842 


15 


. 2688 


. 2942 


.42 


.983 


. 828 


10 


. 3005 


. 3203 


. 41 


.944 


.752 


5 


. 3220 


. 3383 


. 40 


. 842 


. 647 


0 


.3216 


.3381 


.37 


.878 


. 659 


5 


. 3499 


.3614 


.35 


. 829 


.591 


10 


.3711 


.3801 


.37 


. 823 


.639 


15 


.3714 


.3804 


. 39 


. 799 


.548 


20 


. 3752 


. 3815 


.39 


.788 


. 536 


25 


. 3768 


.3817 


.39 


. 642 


.473 


30 


. 3463 


. 3539 


.38 


.582 


.446 


35 


. 3290 


. 3373 


.36 


. 634 


.520 


40 


. 3099 


. 3206 


. 33 


.630 


. 468 


45 


.2914 


. 3011 


. 34 


. 630 


.467 


50 


. 2763 


. 2893 


. 34 


.877 


.674 


55 


. 2665 


. 2793 


. 37 


. 719 


.555 


60 


. 2369 


. 2530 


.31 


. 962 


. 816 


65N 


.1815 


. 2085 


.30 


. 985 


. 796 


Wt avg 


. 3264 


. 3383 


. 3740 


. 796 


. 596 



Table V. Comparison of zonally averaged longwave flux to 

space, FF2, as found by this study for both full- 
CL and 2/3-CL cloud cases for 16 April 1974, and F2(RAS) 
by Raschke et al (1973) based upon NIMBUS III measurements. 
Also included are composite cloud-amount fractions used in 
the two cloud models of this study. Flux values in ly min” 



16 


J anuary 


1974 


16 April 


1974 


16 July 


1974 


Lat F io 


*(f -CL) 


-| -CL 


F io*(f - CL) 


-CL 


F 10*<! -CL) 


f-CL 


20S 


.078 


.758 


.074 


. 842 


. 151 


. 269 


15 


.079 


.733 


.067 


.828 


. 133 


. 306 


10 


.086 


.634 


.071 


. 752 


.135 


. 258 


5 


.090 


. 536 


. 081 


.647 


.093 


.498 


0 


. 082 


. 556 


.073 


. 659 


. 100 


.445 


5 


.086 


.527 


.080 


. 591 


.095 


.534 


10 


.107 


.376 


.086 


. 639 


. 102 


.460 


15 


. 115 


.342 


. 092 


. 548 


. 103 


.430 


20 


.120 


.279 


.085 


.536 


.094 


.431 


25 


. 120 


.320 


. 094 


.473 


.098 


.375 


30 


.111 


. 411 


.098 


.446 


.092 


.425 


35 


. 122 


.382 


. 112 


.520 


.092 


. 383 


40 


. 108 


. 529 


. 116 


.468 


. 085 


.387 


45 


.067 


. 740 


. 100 


.467 


.071 


.566 


50 


.069 


.759 


.080 


.674 


.064 


. 475 


55 


.068 


. 749 


. 108 


.555 


.038 


.819 


60 


. 106 


. 538 


.054 


. 816 


.085 


.459 


65N 


.095 


o 

• 

o 


.082 


. 796 


.061 


. 590 


Wt Avg 


.096 


.514 


.086 


. 596 


.092 


.405 



Table VI. Cross-seasonal comparison of the zonally-averaged 
net fluxes at the surface F^*, and cloud coverage 
based on 2/3-CL parameterization for the six-month interval 

of 16 January 1974 to 16 July 1974. Flux values in ly min" 1 . 
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IV 



SOLAR RADIATION 



A. PARTITION OF SOLAR INSOLATION 

The solar constant assumed in this study at level k=0 
(top of atmosphere) was 2.00 £y rain ^ (Joseph, 1971). This 
flux was depleted by 4% to account for the attenuation 
caused by oxygen and ozone above the tropopause. This left 
the value S = 1.92 £y rain ^ at level k=2 to be used in this 
study as the effective solar constant . 

Eq . (4-1) was then utilized to compute the effective 

solar insolation at the tropopause (k=2) as follows: 

F ( 2 ) = S [^-]“ 2 Cos z (4-1) 

m 

where S - effective solar constant at k=2 

Cos z = cosine of the zenith angle for the Julian date 
used 

r/r = ratio of the actual earth-sun distance to the 

mean earth-sun distance for the Julian date used 
in this study. 

The Smithsonian Meteorological Tables gives the ratio r/r^ 
and the solar declination 6 for 16 April, 0000GMT, as these 
values are employed in Eq. 4-1. These values are, respec- 
tively : 



— = 1.00333 
r 

m 



6 = 8.4833 degrees of lat. 
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6 is used in evaluating the cosine of the solar zenith 
angle, given by 

Cos z = Sin <J> Sin 6 + Cos (p Cos 6 Cos h (4-2) 

where (p is the latitude and h is the hour angle of the sun 
for the meridional data lines. For example. Fig. 1 makes 



that for 




line 1 


h - 55 


line 2 


h = 10 


lines 3,4 


h = 35 



at the times of the synoptic charts (0000GMT and 1200GMT, 
as applicable). Sin <j> was computed according to the standard 
polar stereographic projection formula applicable to the FNWC 
base chart given by 

r 2 - [ (1-32) 2 + ( J-32) 2 ] 

Sin <f> = 5 5 - (4-3) 

r g + [ (1-32 j + (J-32) ^ ] 

2 

where r = 973.752 . Thus Sin <p assumes the following func- 
h 

tional form in terms of the FNWC grid-coordinate I (Fig. 1): 



Lines 1,3,4 



Line 2 



Sin <J> 



Sin (p 



973. 752- 2(32-1) 2 
973.752+2(32-I) 2 

973 . 752- (32-1) 2 
973. 752+(32-I) 2 



(4-4(a)) 



(4-4 (b) ) 
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or, conversely I is given in terms of (p by 



Lines 1,3,4 I = 32 - 22 . 065 [ -■ ] (4-4<c)) 

Line 2 I = 32 - 31. 205 [ j (4-4(d)) 

I = 1,..., 25 for Lines 1,2 

I - 9,..., 25 for Line 3 
I = 63,.#., 38 for Line 4 

For lines 1,2, and 3 the gridpoint soundings correspond to 
0000GMT 16 April when solar noon occurs at the 180th meri- 
dian, Line 4 gridpoint soundings correspond to 1200GMT, 

16 April, when solar noon was at the Greenwich meridian, 

A very simple partition of solar insolation was utilized 
in this study after Joseph (1971), It consisted of dividing 
the insolation F(2) into two parts at level k=2. One part 
was considered to include all wavelengths X > ,9 ym where 
absorption by water vapor and carbon dioxide bands are the 
most prevalent attenuation processes in clear air. This part 
of the solar spectrum was termed the F(A) energy and con- 
sidered subject to water-vapor absorption but not to Rayleigh 
scattering. For those wavelengths A _< • 9 ym, absorption of 

the solar insolation energy by water vapor was considered 
negligible. This part of the solar insolation was denoted 
F(S) suggestive of the fact that it was subject only to Ray- 
leigh scattering attenuation in clear air. The two partitions 
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are formulated after Joseph (1971) as follows: 



F (A) = .349 F ( 2 ) (4-5(a)) 

F(S) = .651 F ( 2 ) . (4-5 (b) ) 

In this study, the introduction of two cloud decks produced 
cloud-reflectivity effects upon both the F(A) and F(S) 
solar energy insolations. However, in the clear areas around 
any gridpoint only the absorption-attenuation applies to the 
F(A) insolation, while only Rayleigh scattering-attenuation 
applies to the F(S) insolation. 

B. DISPOSITION OF F(S) INSOLATION 

In the disposition of the F(S) insolation, Joseph (1971) 
determined that Rayleigh scattering reflectance to space by 
clear skies (after Coulson, 1959) could be effectively approxi- 
mated by least squares in the following form 

a(R) = . 085 + . 25074 [log (j 1 - Sec z)] (4-6) 

o 

where P = 1013.25 mb. In Eq. 4-6, tt/P = 1 in view of the 
o n 9 o 

fact that mean sea level pressure Tf is close to 1000 mb. 

Also 

Sec z = (Cos z) ^ 

with Cos z given by Eq . 4-2. 

The surface albedo a (G) is another reflective parameter 
utilized in this study. Over oceanic areas the following 
formula for a (G) after Gates et al (1971), was utilized? 

a(G) = max {. 06 , . 06 + . 54 (. 7 - Cos z) } . (4-7) 
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1 . Clear Sky Case 



In the clear sky (0,0) case the F(S) insolation was 
subjected to both Rayleigh scattering reflectance a(R) and 
to surface reflectance a (G) • Considering the likelihood of 
a succession of multiple reflections between earth and atmos- 
phere, the F(S) insolation actually penetrating the earth's 
surface after scattering is given by 

IS10(0,0) = F(S) [l-ct(R) ] [l+a(R)a(G)+. . . (a(R)a(G) ) n 

+...]*(l-a(G)) ( 4-8 (a) ) 

that is, by 

IS10(0,0) = F(S) [l-ct(R) ] [l-ot(G) ]/[l-a(R)a(G) ] (4-8(b)) 

2 . Cloudy-Sky Cases 

In the three cases in which clouds were present, F(S) 
insolation absorbed by the ground at each gridpoint was com- 
puted using the following equation (after Arakawa, 1972): 

IS10 ( 1 , 1) = F(S) (l-R(l)) (1-R(2)) (l-a(G)) 

* l-[R(l)R(2)+R(2)a(G)+R(l)a(G) 

+ 2R(l)R(2)a(G) ] -1 . (4-9) 

As indicated by the notation ( 1 , 1 ) , d eno t ing CL ( 1) =CL ( 2 ) = 1 . 0 y 
Equation 4-9 is the formula used in calculating F(S) insola- 
tion absorbed by the ground in the case where overcast clouds 
are present at both levels of Fig. 2. Also in Eq . 4-9, con- 

stant cloud-reflectance values were chosen, namely R(l) = «54 
for the mid-level clouds between k=4 and 6 , and R(2) = 066 
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for the low-level clouds between k=8 and 9. Both cloud- 
reflectance values are as suggested by C. D. Rodgers (1967). 

For all the other cloud cases, the following changes 
were applied to Eq. 4-9. In the (1,0) case (CL(1) = 1.0, 

CL(2) = 0.0), the desired result is obtained by setting 

R(2) = 0 in (4-9), from which it follows that 

I S 10 ( 1 , 0 ) = F(S) (l-R(l) ) (l-a(G) ) / [l-R(l)a(G) ] • (4-10) 

In the case (0,1), one sets R(l) = 0.0 in (4-9) so that (4-9) 
simplifies to 

IS10(0,1) = F(S) (1-R(2) ) (l-a(G) )/ [l-R(2)a(G) ] . (4-11) 

Note that with a cloud overcast present, the Rayleigh clear- 
sky s ca t ter ing a (R) does not appear in Eqs . 4-9, 4-10 or 4-11, 

but is included in the cloud reflectances R(l) and/or R(2). 

3 . Composite F(S) Insolation 

Eqs. 4-8, 4-9, 4-10, and 4-11 were utilized in the 
computation of the cloud-weighted F(S) insolation penetrating 
the earth’s surface considering the areal-weights of the cloud 
combinations denoted by (0,0), (1,1), (1,0) and (0,1) about 

a gridpoint. The resultant F(S) insolation penetrating the 
earth’s surface denoted by IS10 is therefore expressible as 

IS10(CL(1) ,CL(2) ) = IS10(0,0) W(0,0) 

+ IS10(1,1) W(l,l) 

+ IS10(1,0) W(1,0) (4-12) 

+ IS10(0, 1) W(0,1) . 
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Here the weighting factors W(0,0), W(l,l), W(1,0) and W(0,1) 
are computed in Eqs. 2-6a, 2-6b, 2-6c, and 2-6d respectively 
Note finally that the part of F(S) insolation reflected to 
space is found by subtracting IS10(CL(1), CL(2)) from F(S). 
Table VII lists the results of radiative sounding as it 
appears for the disposition of F ( S ) - ins ola t i on . The indivi- 
dual computations of IS10 as they apply for the possible 
overcast-clear layer cases are made under the heading "IS10. 
The difference 

REFS = F ( S ) - IS10 (4-13) 

in each case represents F (S ) -ins ola t i on reflected to space 
wh ile 

STRAK ■ (l E - S a(G)) (4 - 14) 

has been computed as that portion of the F ( S ) - ins o la t ion 
incident at the sea surface just prior to transmission by 
the surface. Note that no absorption in air has been in- 
cluded in the computations of Table VII, and that the only 
absorption permitted is that implicit in IS10. Finally at 
the bottom of each column, e.g., IS10, the composite value 

has been computed by means of the weighting scheme of 
Eq. 4-12. 
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CL ( 1 ) , CL ( 2 ) Weighting 

0,0 0.0 

1,0 0.0 

0,1 .1013 

1,1 .8997 

Composite F(S) values 

Table VII. A sample listing of 
(ly min”-*-) computed 
equations detailed in Sec. IV. 



IS10 


REFS 


STRAN 


. 4305 


.1751 


. 5216 


. 2539 


.3517 


. 3076 


.1921 


.4135 


. 2327 


.1340 


. 4656 


.1696 


.1452 


.4603 


. 1760 



values of F(S) insolation 
at gridpoint (1,1) using 



C. DISPOSITION OF F (A) INSOLATION 

The fractional portion of the solar insolation subject 
to absorption by atmospheric water-vapor and carbon dioxide 
are covered in the following subsections. 

1 • Clear-Sky Case (0,0) 

The Manabe-Moller absorptivity function provided the 
necessary absorptivity values for the key layers in this 
case. The form of this absorptivity function is 

3 0 3 

a(2,k) = .271[U(2,k) Sec z] #JUJ . (4-15) 



Here ji is the absorptivity applied to the pressure-scaled 
water vapor mass between levels 2 and k (Fig. 2) along the 
zenith slant-path angle z. The resultant absorbed insola- 
tional energy in the particular layer (2-6) is then given 
by the Manabe-Moller relation 

30 3 

A ( 2 , 6 ) = 0 . 2 7 IF ( A) [ U ( 2 , 6 ) Sec z]* . (4-16) 
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The two layers of interest in which absorption was computed 
were (2,6) and (2,10). The absorbed insolation in the layer 
(6,10) was then computed by 

A ( 6 , 10) = A ( 2 , 10 ) - A ( 2 , 6 ) . (4-17) 

Water-vapor mass above level 2 was assumed negligible in the 
F(A) disposition of the solar insolation. 

By subtracting A(2,10) from F(A), the direct trans- 
mission of F(A) insolation impinging at the earth's surface 
was determined. The transmission of F(A) insolation is then 
further reduced by the transmissivity (l-a(G)), after sur- 
face-reflectance which leads to the earth-absorbed result 

'JAQ 

IA10 (0,0) = F (A) {l-.271[U(2,10)Sec z ] ° } ( 1-a (G) ) . 

(4-18) 

The transmitted energy impinging upon the earth just prior 
to absorption is 

TRANA (0,0) = IA10(0,0)/[l-a(G) ] . (4-19) 

2 . Cloudy Cases 

In order to compute meaningful dispositions of F(A) 
insolation in cloudy-sky cases, cloud reflectivities and 
cloud absorp tivities after C. D. Rodgers (1967) were utilized. 
The reflectance-values used here are different from those 
suggested by Rodgers for the F(S) wavelengths. The cloud 
reflectivities used here are RA(1) = .46 and RA(2) = .50. 
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In this case there are also cloud-absorp tivi ties to be con- 
sidered. These were taken as A(l) = .20 and A(2) = .30, 
respectively. In the following discussions the cloud con- 
ditions are considered totally overcast in each of layers 
indicated when the notations ((1,1), (0,1) and (1,0)) are 

utilized . 

Schematic representation of the computations per- 
formed in the various cases ((1,1), (0,1), (1,0)) are dis- 

played in Fig. 3. This figure indicates the theoretical 
parameters required in computing the disposition of incoming 
F (A) -ins ola tion from level k=2 to the earth’s surface (k=10) . 
The equations which relate to the parameters in Fig. 3 and 
the other cloud configurations are listed in Appendix B of 
Spaeth’s study (1975). 

The insolations, A24, A46, A68, A89 and A910, etc., 
represent the contributions to the insolation absorbed in 
the layers (2,6) and (6,10) involved. Symbols F2, F4, F6, 

F8 and F9, etc., depict the streams of insolation passing 
through the indicated level. A vertical arrow implies the 
direction if insolation passage, i.e., 1 denotes downward 

insolation, + upwa r d- r e f 1 ec te d insolation, and downward- 

reflected insolation. The absorption quantity A(6,8)ll, for 
example, indicates absorbed energy remaining in (6,8) from 
a downward reflected beam. 

Since multiple reflections occur between the earth’s 
surface and a cloud base, or between cloud-layers in this 
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0 



F(A) 



A24 1 AJ44 





A910U A910T A910* 



10 £ MOT \ / not 

//<»*•/////////////////// 

IA10 

Figure J. Schematic representation of F(A) inso].ation 
disposition in the case of two overcast layers(c a s e(l,1)). 
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model, the insolational amount reaching the lowermost reflect- 
ing surface is affected by the downward reflected beams. 
Additional insolational amounts remaining after more than 
two reflections are negligible, and are not included in the 
model computations. Thus, note that insolation reflected up- 
ward from a lower interface, either cloud-top or ground, to 
the base of a higher cloud deck was not subjected to absorp- 
tion by or transmission through the cloud. This simplifica- 
tion resulted in slightly reduced F (A) insolation-reflectance 
to space, 

A final point to note is that there are three con- 
tributions to the absorption of solar insolation within a 
layer between any two reflecting surfaces. In Fig, 3 these 
contributions are identified by arrows which indicate the 
portion of path being crossed. 

For each of the possible overcast-cloud combinations, 
the transmitted F(A) insolation arriving at the earth's 
surface may be defined using the following notation (TRANA) 
as 



TRANA (1,1) = F104- 


+ 


Ficrn 


(4-20) 


TRANA (1,0) = F104- 


+ 


F104-4- 


(4-21) 


TRANA (0,1) = F104- 


+ 


F1014- 


(4-22) 



where the right side parameters of (4-20), (4-21), and (4-22) 

are derived in Appendix B of Spaeth (1975), 
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The F(A) insolation absorbed by the earth in each 
cloud-case considered is derived from 

IA10 (1,1) = F 10 + ( 1-a (G) ) + FlO't-4- (4-23) 

IA10 (1,0) = F10+ (l-a(G) ) + F10 + 4- (4-24) 

IA 10(0,1) = F104- ( 1-a (G) ) + F104--1 (4-25) 

In the last three equations, the quantity F104-4- is small 
enough in each case, that no further reflections from the 
earth were considered. 

3 . Composite F(A) Layer-Absorptions and Surface - 
Absorption Insolation 

As has been previously discussed, the standard grid- 
area weighting scheme of this study was applied to obtain 
composite values of the absorbed F ( A) -ins ol ati on in key 
layers and also within the earth's surface. The weighting 
factors applied to the corresponding overcast-combination 
absorption quantities provided the following composite results 

A26 (CL ( 1) , CL ( 2 ) ) = A 26 ( 0 , 0 ) W ( 0 , 0 ) +A2 6 ( 1 , 0 ) W ( 1 , 0 ) 

+ A26 (0, 1)W(0, 1) (4-26) 

+ A2 6 ( 1 , 1 ) W ( 1 , 1) 

A6 10 ( CL ( 1 ) , CL ( 2 ) ) =A 610(0, 0)W(0,0) +A 610(1, 0)W (1,0) 

+ A 610(0, 1)W(0,1) (4-27) 

+A 610(1, 1 ) W ( 1 , 1) 

IA10 (CL ( 1) , CL ( 2 ) ) =IA10(0,0)W(0,0)+IA10(1,0)W(1,0) 

+IA10 (0, 1)W(0, 1) (4-28) 

+IA10 (1, 1)W (1, 1) 
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The weighting factors W ( 0 , 0 W ( 1 , 1 ) were first listed in 
Eqs. 2-6 (a,b, c, d) . Also A26(0,0), A610(0,0), IA10(0,0) are 
given in each clear-sky case (0,0) about each gridpoint by 
Eqd. 4-16, 4-17 and 4-18, respectively. 

The results for the absorption in layers (2,6), (6, 

10) and at the surface are shown in the following table for 
the four cloud cases [(0,0), (1,0), (0,1), (1,1)] and for the 
composite cases based on the weighting factors (also listed). 



CL ( 1) , CL ( 2 ) 


Weighting 


A 2 6 


A610 


AI10 


REFA 


TRANA 


0,0 


0.0 


.0809 


.0804 


.1348 


.1751 


. 7069 


1,0 


0.0 


.1192 


.0436 


.0482 


. 3517 


. 37 74 


0,1 


.1013 


.0845 


.1499 


.0238 


.4140 


. 2697 


1,1 


. 8987 


. 1192 


.0717 


.0104 


.4656 


.1880 


Composite values 




. 1157 


.0797 


.0118 


. 4603 


.1962 



Table VIII. A sample listing of gridpoint values (ly (min) ) 
of F(A) insolation computed at gridpoint (1,1) 
using Eqs. 4-26, 4-27 and 4-28. 

In the computational scheme indicated by the entries 
of Table VIII, the reflected F(A) insolation to space has 
been depicted by the symbol REFA, and its values follow from 

REFA = F ( A ) - A 2 6 - A610 - IA10 (4-29) 



whereas the TRANA dispositions are given by Eqs. 4-19,..., 
4-22, respectively, or by its weighted-mean value in the 
case of TRANA-compos i te . 
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4 • Absorptivity (ABA) by Layers 



Here the (fractional) absorptivity as well as the 
actual insolation values absorbed in the layers are con- 
sidered. In the computation of absorptivity, which is frac- 
tional absorption, the total undepleted insolation at the 
top (k=0) is used as a base. The following equation was 
utilized in this calculation: 

FAD J = 2 . 00 (r /r )~ 2 Cos z . (4-30) 

m 



The absorptivity of the troposphere ABA was computed from 
the ratio of the insolation absorbed in the troposphere to 
the insolation incident at the top of the atmosphere rather 
than at k=2: 



ABA 



A ( 2 . 6 ) + A(6. 10) 
FAD J 



(4-31) 



D. ALBEDO (ALB) OF THE EARTH-TROPOSPHERE SYSTEM 

In considering the planetary albedo, the reflected inso- 
lations of the ear th- t r opos phe r e system in both the F(A) and 
F(S) insolational regions must be recalled by the program. 
Thus REF is computed at each gridpoint as the sum of the re- 
flected insolation energy in F (A) , denoted REFA in Eq . 4-29 
and the reflected part of F(S) (denoted REFS in Eq. 4-13): 



REF = REFS + REFA . (4-32) 

Finally the planetary albedo is related to FADJ through 



ALB 



REF 

FADJ * 



(4-33) 
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E. COMPOSITE ABSORPTIVITY ( ABG ) BY THE EARTH-SURFACE; 

COMPOSITE ATMOSPHERIC TRANSMISSIVITY (ATRAN) 

!• Absorptivity (ABG) of Earth 

By summing the weighted values of F(S) and F(A) por- 
tions of the incoming insolation entering the earth, the 
total insolation absorbed at the earth’s surface was com- 
puted. This quantity when divided by the extraterrestrial 
insolation gave the fractional absorptivity (ABG) of the 
earth’s surface. The equation for ABG was 



ABG 



I A 10 + IS10 
FAD J 



(4-34) 



where IA10, IS10, and FADJ were defined previously by Eqs. 

4-28, 4-12, and 4-30, respectively. 

2 . Transmissivity (ATRAN) of the Troposphere 

Also computed was the total insolational energy TRAN, 
incident at the earth’s surface just before absorption by 
the surface. This calculation is given by 

TRAN = TRANA + STRAN . (4-35) 

Here STRAN = [ I S 10/ ( 1-a (G ) ] was previously defined in Eq. 

4-14, and a(G) was given in Eq . 4-7. TRANA has also been 
defined as the weighted value of TRANA(0,0), TRANA (1,1), 
TRANA (1,0) and TRANA (0,1) given by Eqs. 4-19, 4-20, 4-21, 
and 4-22. Note also in justification of STRAN that the four 
cases for IS10 of (4-8), (4-9), (4-10), and (4-11) each have 

the common factor (l-a(G)) in the numerator and therefore 
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each transmitted F(S) insolation component available at the 
earth just before absorption needs only be divided by 
(l-a(G,)). TRAN may thus be viewed as the total insolational 
energy incident at a pyrheliome ter located at earth. The 
( f r ac t ional) tr ansmiss ivi ty of the troposphere (ATRAN) is 
then computed from 



ATRAN = TRAN/ FAD J . (4-36) 

Note finally that the major dispositions of the total in- 
solation at the indicated map times have now been identified 
by the fractional values, ALB, ABA or ABG, and ATRAN, repre- 
senting the reflectivity (albedo), absorptivity of air or 
earth, or atmospheric transmissivity as the case may be, 

3 • Computational Check 

The computational scheme utilized in this model was 
checked continuously by summing the fractional values ALB, 
ABA, and ABG at each gridpoint. The value in each case must 
total .96, since as previously noted the attenuation of 
solar insolation was taken as four percent as it passed 
through the stratosphere. 

F. STATISTICAL ANALYSIS 

In order to substantiate some of the computations per- 
formed in this section, several of the most important items 
computed were statistically analyzed using linear regres- 
sion computer programs from the BMD set of statistical pro- 
grams (Dixon, 1973). 
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1 • Clear Sky Cases 



Using ALB (0,0, z), which is the clear sky case of 
albedo, as the predictand, and log 1Q Sec z as the predictor, 
the following best-fit equation resulted 

ALB (0,0, z ) = . 07734 + . 400721og 10 Sec z 

R m = .9624 (4-37) 

where represents the multiple correlation coefficient. 

This result was as expected over the ocean where 
ot(G) and a(R) involved logarithmic dependence on Sec z. The 
sample average values of ALB and log^ Sec z were 

ALB (0,0) = .1358 log 1 Stec z = .1459 . 

Other clear-sky regression tests made use of the water-vapor 
mass path length (M) defined by 

1/2 

M = (U Sec z) log^^(U Sec z) . (4-38) 

This parameterization of water vapor mass is similar to that 
developed by Hanson (1971), who used a similar M as a pre- 
dictor in his empirical formulations of ABA for both clear 
and partly cloudy sky cases. In this study the two regres- 
sions attempted using M as a predictor were with ABA(0,0,M) 
and ATRAN(0,0,M) as predictands. The best-fit equations 
resulting were 

ABA ( 0 , 0 , M) = .09572 + .03136M (4-39) 

R m = .9875 
M 

AT RAN ( 0 , 0 , M) = .79520 - .03182M (4-40) 

R = .9488 . 

M 
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The absorptivity in the form given by (4-39) agrees closely 
in form with that obtained by Spaeth based on 16 January 
1974 data and with that of Quinn based on pyrheliome tr i c 
data (1971). 

The means of ABA(0,0,M) and ATRAN(0,0,M) for the 
original 93 gridpoint sample were 

ABA ( 0 , 0 , M) = .13459 
ATRAN ( 0 , 0 , M) = .75576 

M = 1.23945 (gm cm" 2 ) 1/2 . 

2 • Statistical Relationships Between ALB, ABA and 

ATRAN in the Cloudy and Clear Sky Cases 

In addition to computing the clear-sky value of ALB, 
it was also possible to compute the composite cloud value of 
ALB. A regression was formed showing the relationship be- 
tween the ratio ALB ( CL ( 1 ) , CL ( 2 ) ) /ALB ( 0 , 0 ) as the predictand 
and the total opaque cloud cover (CL) as the predictor. CL 
in this model is specified by Eq . 3-12. CL, it was felt, 
gave a good approximation of the effective cloud-cover by 
the two layers of cloud amounts CL^ and CL^ by Eq. 2-4 used 
in this model. As explained earlier in this study, a one- 
third reduction of both CL(1) and CL(2) after initial determi- 
nation by Eq . 2-4 was also tried in all computations involving 

cloud amounts. This test represented an attempt at tuning 
the cloud model for radiative calculations, and gave results 
in closer agreement with the latest literature on satellite 
reflectances. Thus all subsequent regressions of ALB, etc. 
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upon CL, after Eq. 2-4, will also include those best-fit 
equations formulated using a CL 1 based upon Eqs, 2-5 and 
3^12. A small letter subscript T a T in the equation number 
will indicate those computations using 2/3-CL amounts while 
a small letter subscript f b f will indicate those utilizing 
full-CL cloud model amounts. 

The first regression equations tested were those for 
ALB(CL) and ABA (CL) for the respective cases of 2/3-CL and 
full-CL cases. The results are given in (4-41( a ))> (4 — 41(b)) • 

The symbol once again signifies the multiple correlation 
of the statistical regressions. 



ALB ( CL ( 1 ) CL ( 2 ) ) = ALB(0,0)[1 + 4.7154CL 
a 

- 2.3648CL 2 ] (4-41 (a) ) 

R m = .9640 
M 

ALB, (CL(1)CL(2)) = ALB(0,0)[1 + 4.7273CL 
b 

- 1.6974CL 2 ] (4-41 (b) ) 

R m = .9670 
M 

The sample means of the values in Eqs. 4-41a and 4-41b are 



Eq . 4-41a 

ALB (CL (1) , CL (2) ) = 



ALB (0,0) 

cT 



. 3624 
.1358 
. 5876 



Eq . 4-4 lb 

ALB (CL (1) , CL ( 2 ) ) = 



ALB (0,0) 

cT 



. 4473 
.1358 
. 7746 



The high value of the correlation coefficient in Eq. 4-41 
shows strong dependence of ALB upon CL in either cloud 
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parameterization case. Also note that ALB ( CL ( 1 ) , CL ( 2 ) ) 
values computed here are clearly higher than reported by 
Raschke et al, (1973), where ALB = .25 for essentially the 
same set of gridpoints during the NIMBUS III period 1-15 
May 1969. The 2/3-CL cloud parameterization gives albedo 
results which are closer to Raschke* s results, suggesting 
that the full-CL amounts of Smagorinsky (Eq. 2-4) are too 
high for radiational computations. 

A second regression was developed between ABA(CL) 
/ABA(0,0) versus CL. The following best-fit equations 
results: 

ABA (CL ( 1) , CL ( 2 ) ) = ABA(0,0)[1 + . 5873CL 
a 

- .1018CL 2 ] (4-42(a)) 

R m = .8956 

ABA,. (CL (1) , CL (2) ) = ABA(0,0)[1 + . 6526CL 
b 

- . 1481CL 2 ] (4 — 42(b)) 

R = .9012 
M 

Mean statistics for this regression case are 
Eq. 4-42a Eq . 4-42b 

ABA (CL (1) ,CL (2) ) = . 1758 ABA, (CL (1) ,CL (2) = .1894 
a b 

ABA (0,0) = .1346 ABA (0,0) = .1346 

CL =. 5876 CL = . 7746 

a b 

From the means of (4-42), it is seen that the model specifies 
an increase in atmospheric (solar) absorptivity with increas- 
ing cloud cover. This result is in agreement with studies 
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of Warner (1974) and Spaeth (1975) and with the earlier 
study of London (1957). 

An analogous statistical regression was then de- 
veloped for the cloudy-sky transmissivity relative to the 
clear-sky transmissivity. The resulting regressions were 

ATRAN (CL (1) , CL ( 2 ) ) = ATRAN(0,0)[1 - . 6581CL] 

R m = . 9992 . (4-43 (a) ) 

ATRAN fe (CL (1) ,CL(2) ) = ATRAN(0,0)[1 - .6833CL] 

R m = .9986 (4-43 (b)) 

with sample mean values as follows: 

Eq. 4-43a Eq . 4-43b 

ATRAN (CL(1) , CL ( 2 ) ) = .4653 ATRAN (CL(1) ,CL(2))= . 3589 
a b 

ATRAN (0, 0) = . 7558 ATRAN (0,0) = . 7558 

The parameter ATRAN (CL ’) /ATRAN (0 , 0 ) given by (4-43a) is in 
very good agreement with the statistical result of Savino- 
Angstrom from earth-based pyrheliometry at Canton Island, 
which was 



Q = Q q (1 - . 6 5 5CL ) 

(Budyko, 1958). The coefficient of CL in this result is 
known to be a slowly decreasing function of latitude. 

Similarly, the parameters ALB (CL ' ) / ALB ( 0 , 0 ) and 
ABA (CL ' ) / ABA ( 0 , 0 ) behave consistently with respect to increas- 
ing CL, with more realistic mean values for ALB(CL') and 
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ABA(CL') chan for the corresponding full-CL parameterization. 
Hence, in the mean, the statistical analyses performed seem 
to show agreement with observational results of other investi- 
gators. Comparative computations of the atmospheric trans- 
missivity, particularly in the 2/3-CL case appear to give 
better agreement with ground-based pyrheliometry . 

A major test still to be made is the comparison of 
the latitudinal-distribution of ALB (CL) and ALB (CL) with 

3 . D 

that of the satellite albedo of Raschke et al (1973). 

G. ALBEDO COMPARISONS WITH PUBLISHED RESULTS 

The tropospheric albedo computations of the solar-inso- 
lation model of this section have been presented for both 
the full-cloud and two-thirds CL cases, respectively. These 
albedos have been interpolated to whole multiples of 5° lati- 
tude between 20°S and 65°N. These computations were made 
for each of the four oceanic meridians and the resulting 
albedos averaged across the meridians are presented as a 
function of latitude in Table IX for both the full-CL and 
two-thirds CL cases, where comparison is also made with 
satellite albedos of Raschke et al (1973), for the NIMBUS 
III period 1-15 May 1969. For sake of consistency the 
Raschke albedo data have been extracted from atlas presenta- 
tions to the 5 degree latitude grids, on the same four 
meridians, before applying the same averaging process. 

The model results under the headings "FULL-CL 11 and "2/3- 
CL" correspond to the composite cloud amounts previously 
listed in Table V according to the Eq . 3-12. 
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The most unusual results to be observed in Table IX are 
the small albedo-values (~ .2) between 20S to 20N, after 
Raschke,. as contrasted with the corresponding values (in the 
range *3 to .47) by the 2/3-CL model, and more markedly with 
those by the full-CL model (values in the range .43 to .59). 
The cosine-weighted mean albedos (see Eq . 7-2) are as 

follows : 



FULL-CL, ALB = .46 
2/3-CL, ALB = .37 
RASCHKE-ALB = .25 

At every latitude ((}) <_ 60) under comparison the 2/3-CL 
albedo computation is closer to Raschke’s observations than 
is the full-CL value so that there is no question in accept- 
ing the 2/3-CL over the full-CL case. At latitude <J) = 65, 
there is evidence of added surface-ice reflections contribut- 
ing to the Raschke value of planetary albedo. For example 
at latitude 65N Raschke’s results indicate ALB = .5, whereas 
the 2/3-CL model computations with a rather high cloud cover 
(CL 1 = .79) is .47. Part of this discrepancy with the 
solar radiation model at high latitudes lies in the non- 
inclusion of a surface-ice reflectance parameter. 

To summarize, apart from the region covered by sea-ice, 
there is a marked preference for the 2/3-CL over the full- 
CL model in giving albedos closer to those observed by 
satellite for the same general period (1-15 May 1969) as 
reported by Raschke et al (1973) and Von der Haar and Hanson 
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ALBEDO MODEL 

Cloud-Amount 



La t . 


(Full-CL) (2/3-CL) 


(Ras chke ) 


(Full- 


CL) (2/3- 


20S 


. 568 


.473 


.17 


.985 


.842 


15 


.559 


.462 


.17 


.983 


.828 


10 


.512 


. 418 


.18 


.944 


.752 


5 


.446 


.363 


.20 


.842 


. 647 


0 


.475 


.374 


. 21 


.878 


.659 


5 


.443 


.347 


.22 


.829 


.591 


10 


.438 


.337 


.20 


.823 


.639 


15 


.434 


.333 


.19 


. 799 


.548 


20 


.430 


.331 


.18 


.788 


. 536 


25 


.381 


.314 


.20 


. 642 


.473 


30 


. 351 


. 298 


.24 


.582 


.446 


35 


.361 


.316 


.29 


. 634 


. 520 


40 


.372 


.311 


.28 


.630 


.468 


45 


.414 


.336 


.30 


. 630 


.467 


50 


.499 


.415 


.33 


.877 


.674 


55 


.443 


.383 


.40 


.719 


. 555 


60 


.562 


.493 


.42 


.962 


.816 


6 5N 


.609 


.474 


. 50 


.995 


. 796 


Wt Avg 


. 458 


.373 


.24 


.796 


.596 



Table IX. Comparison of planetary albedo as found by this 
study for both full- and two-thirds cloud models 
for 16 April 1974, and by Raschke et al (1973) based upon 
NIMBUS III measurements. Also included are the globally- 
weighted mean values at the bottom of each tabular column 
and the composite cloud amounts for both cloud cases. 
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(1969) and others. The primary difficulty even with the 
2/3-model occurs in tropical and subtropical latitudes 
where Von der Haar and Hanson (1969) have discussed the 
reality of even smaller than suspected cloud-covers and 
smaller values of the resultant global-albedo. At these 
latitudes (tropical and subtropical) it seems doubtful that 
a large-scale CL-parameteri zation of the Smagorinsky type 
will be applicable. Rather some method of specifying the 
distribution at gridpoints of isolated cumulus towers may 
be appropriate. Finally, it is open to question whether 
vertically-structured convective cloud elements will have 
the high reflective capability attributed to clouds which 
are depicted as existing in horizontal layers with a con- 
stant reflectivity . 5 R <. .65 as specified in the radia- 
tion model. It may well be that the vertical shafts between 
convective elements may be a more efficient focus for solar 
energy to the earth than is perceived by ascribing to the 
cloud amount a constant cloud albedo value. A more realistic 
approach may be to reduce the two-cloud layering in the 
tropics to clouds existing in a single layer. 

At any rate with the two layered cloud-model considered 
in this study, the only feasible choice is the 2/3-CL para- 
meterization. This selection is carried over into the com- 
putations of Section VI. A further test of the 
c 1 oud -pa r ame t er i za t i o n effect on the computation of net 
radiation of the earth-atmosphere system, by comparison with 
that observed by satellite, is made in Section VII. D. 
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V. 



SENSIBLE AND LATENT HEAT TRANSPORT 
AT THE SEA-AIR INTERFACE 



A. GENERAL PURPOSE 

The model used here to describe the turbulent transports 
of sensible and latent heat across the sea-air interface was 
basically that already in use in the operational FNWC primi- 
tive equation model. This model has been discussed in detail 
by Kesel and Winninghoff (1972) and by Kaitala (1974). The 
primary purpose of the adaptation of the turbulent flux model 
in this study was to test the magnitudes and directions of 
these transfer rates in comparison with those of the radia- 
tional transfer model discussed here in Sections III and IV. 
The final purpose of the inclusion of the turbulent flux 
model was to conduct a diagnostic computation of the heat 
balance within the oceanic and atmospheric layers for 16 April 
1974 under the combined operation of the radiational and 
turbulent flux models. 

B. WINDSPEED COMPUTATION IN THE TURBULENT FLUX MODEL 

At all of the 93 gridpoints tested, the geostrophic wind- 
speed Vg was computed by use of the standard formula 




(5-1) 



wi th 




(5-2) 



1.86603 
1 + S in <J) 
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and 



AZ =-y/(Z 4 -Z 2 ) 2 + (Z 3 -Z x ) 2 (5-3) 

In Eq. 5-1, d = 381.0km is the nominal FNWC grid-spacing 

o' 2d 

(true at 60 N) and (-^-) is the true earth spacing used in 
the centered-dif f erence computations of the contour gradient 
( cf . , Fig . 3). 

In Eqs. 5-1 and 5-3 AZ is the contour-difference, cen- 
tered in most cases on the gridpoint under consideration. 
However, at the map edge gridpoints of the meridians, it was 
only possible to compute a f o rwa r d- d i f f e r enc ed Version of Vg 
(cf.. Fig. 4). Thus at edge-point "0", Eq. 5-4 below was 
used. The spacing — in Eq. 5-4 is then taken as one-half 
that in (5-1), so that 

Vg = -2- — A2 (5-4) 

Vg f (d/m) 

with 

AZ = y (Z 4 -Z 0 ) 2 + (Z 3 -Z 0 ) 2 (5-5) 

All geostrophic windspeeds were then converted to surface 
windspeeds using the following empiricism due to Langlois and 
Kwok, (1969) 

V = .8 Vg + 2.2 (5-6) 

s 

where V = surface windspeed (mps) and Vg is the lOOOmb geos- 
s 

trophic windspeed. The factor 2.2 in (5-6) is an empiricism 
introduced to include the effects of gustiness upon the 
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measured (steady-state) windspeed. 

In the interval of + 10 degrees of latitude, the value 
of the Coriolis parameter was arbitrarily set at f =+,25 x 
10 ^sec following Kaitala (19 74), while for | <J> | > 10 de- 
grees the actual value of f for the particular latitude was 
used . 



C. SENSIBLE HEAT TRANSPORT 

The details of the adaptation of the FNWC model for sen- 
sible heat transport Hp utilized in this study follow the 
description of Spaeth (1975)* The formulation of Hp is 
given below in terms of bulk parameters of the surface and 
of the planetary boundary layer (10,8) of Figs. 2 and 5: 



H r " p io C P c d V T g " V 



10 



P 10 (R d T 10 ) 



-1 



(5-7) 



where 



and 



C = .239 cal gm 1 (deg.K) 
P 



C = 1.4 x 10 (after Weiler and Burling, 
D 1967) 



The equation utilized in this study to compute may be 
shown to be (Spaeth, 1975) 



T 

x 



CV T 
D s g 



2 k e 9 

AZ 



2 K 
AZ 



+ cv 

D s 



2 K Y 

c 



(5-8) 
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where 



fl _ T ,1000. .286 
y ‘ T 9 < ToF ) 



K = 



K 



. e -e 

l+a*(=2_S_) 



V Z 10 



AZ = 



Z - Z 
9 10 



f 1 Un 



(5-9) 

(5-10) 



(5-11) 



In Eq. 5-10, the turbulent parameters K* and a* have been 
determined empirically by FNWC 

K* = 10 cm sec 

a* = 5 x 10 4 cm(deg.k) ^ . 

T - T in Eq. 5-7 is the temperature change for the thin 
8 x 

layer between the surface and the top of the constant flux 
layer (cf . , Fig. 5(a)). The formula for deduced here re- 
quired a constant value of sensible heat Hp then decreasing 
linearly with pressure to Hp = 0 at k=8, and remaining zero 
above k=8. This linear decrease of Hp with pressure is com- 
patible with the condition of constant convergence of the 
sensible heat flux per gram in the turbulent boundary layer 
(x,8). Thus, the constant value of the convergence of sen- 
sible heat as given by Eq. 5-7 within the layer (8,10) may 
then be considered to be applicable at level k=9 as depicted 
in Fig. 5(a), after Spaeth (1975). 

After computing T^ by Eq. 5-8, Hp was then computed at 
each gridpoint using the bulk transfer equation 5-7. The 
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Figure k. Section of FNWC polar stereographic meridian il- 
lustrating the method of obtaining contour gradients in the 
vicinity of a gridpoint. The values of contour-heights 
correspond to the 1000 mb D-values in the neighborhood of 
the gridpoint. 
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gridpoint values of Hp are used later in calculations of 
atmospheric and oceanic heat balances in Sections VI and 
VIII. 

D. EVAPORATIVE HEAT TRANSPORT 

Similar to the case of sensible heat convergence, the 
latent heat transport by turbulence was assumed to be sub- 
ject to cons tan t evaporative f lux convergence throughout 
the layer 800-1000mb. The total amount of latent heat re- 
moved by evaporation from the ocean surface was essentially 
as modeled by Kaitala (1974) using b ulk- t r ans f er theory. It 
should be recognized however, that the condensation and con- 
sequent realized latent heat may occur at arbitrary levels, 
i.e,, in general at levels higher than k=8. Thus, in fol- 
lowing the FNWC model, we have not introduced E as a heating 
rate at level k=9. 

The basis of the turbulent latent heat model adopted is 
shown in Fig, 5(c), where 

E (900) = 1/2 E (Bulk) (5-12) 

and where the surface layer bulk transfer of latent heat is 
given by 

E = E ( Bulk) = L P 10 c D v s (<i s " q x > • (5-13) 



Evaporation at level k=9 may also be computed using 

~ (q x " V 

E (900) " LP 10 K (Z 9 -Z 1Q ) 
where K is given in Eq . 5-10. 



(5-14) 
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In (5-13) and (5-14) L is the latent heat of vaporization. 



L = 596.73 - 601T 1q 

in cal(gm) ^ . In (5-13) and (5-14) q and q are the mix- 

9 x 

ing ratios at the levels k=9 and at x, respectively. The 
factor 1/2 which appears in the right side of (5-12) re- 
emphasizes the concept of constant latent-heat convergence 
in the level (x,8), (cf.. Fig. 5(c)). 

The solution for q^ is obtained by combining Eqs. 5-12, 
5-13 and 5-14 in the form 



q 



x 



C D V s q l0 + 2 AZ q 9 




2 



K 

AZ 



(5-15) 



where AZ = - Z . 

Applications of Eqs. 5-13 and 5-15 make possible the 

computation of E at each gridpoint. As in the case for T^, 

the factor of 2K of (5-15) was taken to conform to the value 

K as used by Kaitala (1974) so that the turbulent transport 

results obtained here will be compatible with those computed 

operationally by the FNWC primitive equation model. As 

already noted, this was accomplished in this study by using 

5 2 -1 

a value of 2K* = 2 x 10 cm sec as equivalent to the Re- 
value of Kaitala (1974). Thus computations of E by Eqs. 

5-13 and 5-15 were equivalent to that using Kaitala ’s scheme 
(1974), the only change being that the vertical distribution 
of E implicit in Fig. 5(b) is slightly different from that 
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of the constancy of E everywhere in the planetary boundary 
layer, as postulated by Kaitala (1974), 

E. TURBULENT HEAT TRANSPORTS OVER AN ICE-COVERED OCEAN 
Though this study did not have any data taken from an 

ice-covered ocean, it should be noted that Spaeth (1975) de- 
vised an acceptable method of handling this realistic con- 
sideration, which would be applicable if our computations 
extended beyond the poleward limit of this study which was 
established at <j> = 65N. 

F. LARGE-SCALE TURBULENT HEAT FLUX ACROSS THE SEA-AIR 
INTERFACE 

As shown in Section VI, the combination of E + Hp turns 
out to be a heat loss mechanism at the ocean surface, as 
computed at each gridpoint. The values of Hp which repre- 
sent the sensible heating rates for the atmospheric layer 
(8,10) were computed for each gridpoint along the four 
meridians. The values of Hp in general turned out to be 
small and negative in low latitudes, reflecting the computed 

negative values of T - T , as also reported by Spaeth (1975). 

8 x 

In the mid-latitudes, where T - T > 0, H r was positive and 

g x r 

therefore contributes positively to the net atmospheric tem- 
perature changes in the layer (8,10). The latent heat term 
(E) was dominantly large and positive relative to lip through- 
out this study. Thus, regardless of possible inconsistencies 
in the sign of Hp, combined values of E + Hp seemed to give 
reasonable heat-loss rates at the ocean surface at all lati- 
tudes, as also found by Spaeth (1975). 
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flux after Kaitala (1974) (c). Note 'x' is the level of the top of the 
flux layer. 



In the large-scale sense E + Hp is included as a positive 
contribution to the heat budget of the tropospheric column, 
even though it is not known precisely where the latent heat 
source E will be realized. The use of the equations of con- 
servation of mass, momentum and of water-vapor mass, would 
resolve the latter question, but this was considered beyond 
the scope of the present work. 

The turbulent flux computations for Hp and E outlined in 
the foregoing subsection V # (C). and (D) represent minor modifi- 
cations of the FNWC turbulent boundary layer model which has 
been appended to our radiation model (for the 2 /'3-CL case) 
for heat balance computational consistency tests, that is, 
to determine whether the net radiative flux at the surface is 
of the proper magnitude, relative to E + Hp . 
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VI 



MERIDIONAL CROSS-SECTIONAL DEPICTION 
OF THE HEAT-BALANCE COMPUTATIONS 



A • GENERAL 

The general design of this section is to utilize all of 
the computational concepts discussed in Sections III, IV and 
V in the computations for a single time-step in the heating 
model developed for use in the FNWC prediction model. After 
testing F 2 * and ALB with the corresponding values of Raschke 
et al (1973), it was decided that only the computations by 
the 2/3-CL parameterization would be displayed in the meri- 
dional cross-sections (Figs, 6, 7, 8, 9 and 10), The appropriate 
calculations were performed at each gridpoint of the four 
meridians used for presentation pruposes in this section. 
However, the same type of computations using the full-CL 
parameterization were also made, but not presented in this 
section, 

B. GEOGRAPHICAL REPRESENTATION OF THE HEAT-BALANCE 

DISTRIBUTION 

The FNWC gridpoint processed analyses for 0000GMT, 16 

April 1974 were used at the three Pacific cross-sections, 
while that for 1200GMT, 16 April* 1974, were used for the 

single Atlantic meridian. This was done so that the set of 
gridpoints was considered to be subject to the actual radia- 
tive-transfer calculations involved for these specific 
times. Figure 6 depicts in symbolic language the key to 
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the computational entries in Fig. 7, 8, 9 and 10. This 
symbolic list presents the computations made at each radia- 
tive sounding gridpoint (I,J) having data in the form of 
Table 1(b). The computations proceed from the top of the 
troposphere to the ocean surface and include the sensible 
and latent heat transfers, by the means of the computations 
of Section V. For purposes of climatological data compari- 
son, Figs. 7, 8, 9 and 10 were developed by interpolating 
gridpoint results to integral multiples of 5-degree incre- 
ments. The interpolation routine to this gridpoint spacing 
made use of the Lagrangian cubic interpolation scheme (after 
Spaeth, 1975) . 



Q (I ) 



n (1-2) (1-3) (1-4) (1-1) (1-3) (1-4) 

^1 (1-2) (1-3) (1-4) ^2 (2-1) (2-3) (2-4) 



(6-1) 

, n (I-D (1-2) (1-4) (1-1) (1-2) (1-3) 

^3 (3-1) (3-2) (3-4) ^4 (4-1) (4-2 ) (4-3) 



Finally for ease in reconciling the magnitudes of all heat- 
transfer rates, the time-dependent solar disposition rates 
have been averaged to 24-hourly rates. 



C. EXPLANATION OF SYMBOLIC TERMS 

1 . Cross-Section at Level k=2 (Fig. 6) 

The discussion of all insolation parameters discussed 
previously in Section IV dealt with the specific time of day 
that corresponded to the hour angle h for the instantaneous 
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time t under consideration. The incident solar insolation 
dealt with is then 

F(S) = S ~ 2 Cos z (6-2) 

m 

In order to avoid reference to specific map times t, the in- 
stantaneous solar hour-angles were h = 35, 10, 55, and 35, 

respectively for cross-sections 1, 2, 3 and 4, as depicted 

in Figs. 7, 8, 9 and 10. 

QAVE represents the 24-hour average of F(2) and 
appears as the first input symbol in Fig. 6. Its value is 
considered to be more representative climatologically for 
the data day under consideration than F(2). 

QAVE is derived by the formula 

QAVE . F (2 ) §2§-f (6-3) 

where 

Cos z= [H Sin(J)Sin 6 + Co s(j> Co s 6 S inH ] / tt (6-4) 

H = arcCos [-Tan(j)Tan(S] (6-5) 

Here 6 is the solar declination for 16 April 1974, and H is 
the hour angle at local sunset at latitude <J> . Cos z in (6-3) 
is equal to the 24-hour average cosine of the zenith angle, 

Eq. 4-2. The 24-hour time averaging period yielding QAVE 
gives heating results consistent in magnitude with the terres- 
trial flux divergences, which change only slightly with the 
time of day. The conversion to expected daily averaged 
solar disposition quantities is compatible with the determi- 
nation of a heat budget for the given date (16 April 1974). 
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Other parameters needed for level k=2 are 



QREF = REF(t) (6-6) 

where 

REF ( t ) = F(2)- A26- A610-(IA10+ IS10) (6-7) 

REF(t) is the instantaneous solar reflected insolation at a 
gridpoint and QREF is its 24-hour average, assuming that 
planetary albedo remains constant for the 24-hour period, 
although it was calculated for the indicated solar time t. 
This assumption requires that the cloud amounts computed at 
the indicated synoptic times are representative of the 
entire day. 

The same principle will be used with regard to all 
other solar parameters in the conversion from time-dependent 
values at solar time t to 24-hour averaged values. Superior 
bars ( ) are not used in presenting the averaged values 

shown in the cross-sections key, Fig. 6. The 24-hour aver- 
age system balance, BALT, is computed from 

BALT = QAVE - (QREF + F 2 *) (6-8) 

for the tropopause level k=2 at the indicated latitude. Net 
terrestrial fluxes, such as F 2 *, were considered to be con- 
stant throughout the 24-hour period, a valid assumption if 
the cloud cover remains quasi-constant for the period. 
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2 



Cross-Section in Layer (2,6) 



The following definitions apply for the layer (2,6) 
as referred to in Fig. 6. All the heat transfers shown in 
this layer are assumed to be of radiative character only, as 
in the corresponding version of the FNWC heating model. The 
averaged radiative heating (cooling) rate is given by 

BAL26 = Q26 - F26 (6-9) 

where 

Q26 = daily solar absorption in layer (2,6) 

and is defined relative to A26(t) by a cosine transformation 
similar to Eq. 6-3, and 

F26 = terrestrial cooling rate. 

3 . Cross-Section in Layer (6,10) 

The 24-hour average radiative cooling in the layer 

(6.10) is given by 

BAL 610 = Q610 - F610 (6-10) 

where Q610 is the sum of Q68 and Q810. BAL 68 has been taken 
as one-half of BAL610. However, the heat balance in (8,10) 
has been modified for inclusion of Hp, as given by Eq. 5-7. 

BAL810 = Q 8 10 - F810 + Hp = BAL 6 8 + Hp (6-11) 

Here Hp is the sensible heat flux convergence in the layer 

(8.10) . F 6 1 0 is the sum of F68 and F810, the terrestrial 
flux losses for the respective layers. 
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4 • Cross-Section at Air-Sea Interface (k=10) 

The heat balance at the earth's surface (BALB) is as 
defined in the following equation 

BALB = QABG - F - (H f + E) (6-12) 

QAGB is the 24-hour average insolation absorbed by the sur- 
face as 

QABG = QABG(t) (Cos" z/Cos z) (6-13) 

Hp and E were computed respectively from Eqs. 5-7 and 5-13 
and considered constant for the computational day, 16 April 
1974 . 

D. LATITUDINAL CROSS-SECTIONS OF THE VERTICAL HEAT 
BALANCE FOR 16 APRIL 1974 

Figs. 7, 8, 9 and 10 as previously explained represent 

the single time step of heating computations for each of the 
four meridians used in this study. The four figures have 
been divided into (a) tropical results and (b) mid-to-high- 
latitude results. While the results depicted in these cross- 
sections are exhibited as representing daily-averaged values, 
they are actually based upon heat-computations at the speci- 
fic map times of 0000GMT and 1200GMT on 16 April 1974. There- 
fore, for these results to be meaningful as a stepwise part 
of the heat package subroutine of FNWC, the solar radiative 
absorption and reflectance terms would have to be 
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recoverable as a function of GMT, i.e<>, 



F ( 2 , t ) = QAVE* (Cos z/Cos z) (6-14) 

REF(t) = QREF* (Cos z/Cos z) (6-15) 

etc. Thus solar disposition terms may then be utilized in 
connection with the one-hour stepwise application of the 
thermodynamic equation of the set of primitive equations 
used in the FNWC prediction process, assuming the 2/3-CL 
parameterization. 

The full-CL heat-balance cross-sections were also com- 
puted, but have not been presented here for the sake of 
brevity. However, certain key radiative heating (cooling) 
rates have been summarized in zonally-averaged form from 
both cloud parameter izations for the purposes of presenting 
the mean latitudinal comparisons (Section VII) of the 
purely radiative heat budget of the FNWC data as it was 
applicable to 16 April 1974. 
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The reduced cloud amounts 2/3 CL(1) and CL(2) are given by Eq 
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for key. Values computed from data for 16 april 197^+ for 2/3-CL case only. 
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Figure 8(a). 170W longitudinal cross-section, tropical section. Refer to Fig. 

for key. Values computed from data for 16 April 197 ^» 
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for key. Values computed from data for 16 April 197 ^* 
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for key. Values computed from data for 16 April 197^* 



5923 0.5752 0. 5540 0.5288 0.4999 0.4675 
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Fig. 6 for key. Values computed from data for 16 April 197^- • 
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for key. Values computed from data for 16 April 197^« 
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VII . THE LATITUDINAL DISTRIBUTION OF 
RADI ATIONAL BALANCE TERMS OF THE 
OCEAN-ATMOSPHERE SYSTEM 



A. GENERAL 

The la t i t ud ina 1 ly-d i s t r ib u t ed results of the purely 
radiative contributions, over the ocean-atmosphere system, 
are presented in Figs* 12(a,b) for the 2/3-CL and 13(a,b) 
for the full-CL cases respectively. These cross-sections 
show the results after averaging over * the four meridians 
considered in this study. The results are displayed in the 
format of Fig. 11. Radiative transfer Q68 and Q810 have 
been combined with F68 and F810 of Fig. 6 to compute layer- 
average radiational warming rates in the layer (6,10). The 
turbulent transfer heating rates E + Hp into the atmosphere 
(and the corresponding cooling rates in the ocean) which 
were included in the results of Figs. 7, 8, 9 and 10 are now 

omitted in the radiative composites presented in Figs. 12 
and 13. 

In obtaining la ti tudinally-dis tr ibut ed means of radia- 
tive heating rates at each five-degree multiple of latitude 
<p > all values of each radiative parameter listed in Fig. 11 
at latitude <J> in the range 20S,...,65N (5-degree increments) 
were simply averaged over the four meridians providing the 
averaged value at <f> • At <p = 65N there was only one contri- 
bution to the latitudinal average, while in the Southern 
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Hemisphere latitudes, 20S, 15 and 10S, only two values (on 

X = 125W and X = 35W) of each parameter contributed to the 
means listed in Figs. 12 and 13. At 5S, there were three 
sets of radiative parameters entering into Q(<j>). Otherwise, 
there were four values of each radiative parameter entering 
into the computed radiative mean values Q(<J>) of the cross- 
sections. Therefore, near the northern and southern 
boundaries of Figs. 12 and 13, the listed values may not be 
equally representative, although the general trend should be 
reliable. Fig. 14, using data extracted from Figs. 12 and 
13, tends to confirm this statement. 

Figs. 12, 13 and 14 are also useful in determining the 

relative merits of the 2/3-CL model as compared to the full- 
CL model. 

B. EARTH-TROPOSPHERE SYSTEM RADI ATI ONAL BALANCE SUMMARY 

The latitudinally-averaged distributions of the three 

key radiative parameters, R , R and R as functions of lati- 

S 3. 

tude for both the 2/3-CL and full-CL cases are shown in 

Fig. 14. The radiative-balance parameters were defined 

(after Malkus , 1962) as R ((j>) , the mean radiative-energy 

s 

transfer rate across the top of the ocean-tropospheric 

system at k=2 (referred to as BALT in Fig. 11); R (0), the 

a 

mean radiative cooling rate in the troposphere (BAL26 + 
BAL610 in Fig. 11); and R(cj>), the mean radiative warming 
(cooling) rate at the earth’s surface (BAL10 in Fig. 11). 
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LAT latitude of the gridpoints in 5 increments 

a) QAVE 24-hour averaged mean insolation at level k=2 

b) QREF mean reflected average insolation at level k=2 

c) F * mean net outgoing long-wave flux at level k=2 

d) BALT mean averaged earth tropospheric gain rate (a-b-c) 
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Figure 12(a). Zonally-averaged radiational cross-section for tropical latitudes with 
2/3-CL model. Refer to Fig. 11 for key. All values listed are daily averages in 
ly min and are computed from data for 16 April 197^* 
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Figure 12(b). Zonally-averaged radiational cross-section for higher latitudes 
with 2/3-CL model. Refer to Fig. 11 for key. All values listed are daily averages 
in ly min -1 and are from data for 16 April 1974. 
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Figure 13(a)* Similar to Fig. 12(a) except full-CL model (after Smagorinsky 
is used in all radiative calculations. 
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igure 13(b). Similar to Fig. 12(b) except full-CL model (after Smagorinsky, 
q. 2-4) is used in all radiative calculations. 



The relationship between these radiative parameters is 



R = R + R 
s a 



(7-1) 



Tables X, XI and XII list the zonally-averaged values of R , 

s 

R and R as a function of latitude <J> for the 2/3-CL and the 

a 

full-CL cases. Weighted averages for any Q-value are then 
computed using the following cosine weighting scheme (Eq. 7-2, 
Spaeth, 1975). 



wt. avg, 



18 / 4 

, i=l v i=l 1 i 



18 



E -t 



Cos <f> 



i=l 



(7-2) 



where Q.. = value on meridian i at latitude d> . and k. is the 

number k=l,...,4 of meridional observations available for the 

meridional averages at <p . Note that i=l,...,18 corresponds 

to the 18 latitudes, <J>^ = 20,. .., 65 latitude degrees. Tables 

X, XI and XII depict comparative values of R , R and R by 

s a 

latitudes for both 2/3-CL and full-CL parame terizat ions for 
16 January, 16 April and 16 July 1974. 

In Fig. 14, the dominant net-flux terms across inter- 
faces are R ((p) and R ( cf) ) . The net flux values obtained from 
s 

the 2/3-CL model (denoted R -) exceeded those by the full-CL 

s 1 

model (R _ ) at all latitudes. Both R functions had maxima 
s 2 s 

at 30N on 16 April 1974. Comparative values of the mean 
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latitudinal distributions of the radiative balance R ( cf) ) at 
sea level for the 2/3-CL and full-CL cases were also ex- 
tracted from Figs. 12 and 13. They were tabulated in the 
center portion of Tables X, XI and XII for 16 April 1974 and 
plotted also as a function of latitude in Fig. 14. It is to 
be noted that the 2/3-CL case (denoted R^) is greater than 
over the entire range of latitude under consideration in 
Fig. 14. The difference R^ - R g ^ is very nearly equal to 
that of R^ - R s 2’ S ° t ^ iat t * ie net cooling in air, R = R £ 
over the entire range of latitude in Fig. 14. 

In summary, for 16 April 1974, Table X and Fig. 14 show 
that a substantial positive difference in net flux exists at 
the troposphere k=2 for the case of R g ^ (for 2/3-CL) over 
the case of R £ (for full-CL). The zonally-aver aged R g -dif- 

ference for the range of latitudes in Spring is .0380 cal 

-2 -1 

cm min • Likewise, Table XI and Fig. 14 show that the net 

flux at the surface for the 2/3-CL exceeds the full-CL by 

-2 -1 

.0351 cal cm min , when both sets of values R^ and R^ have 
been meridionally-averaged. However, the difference ^ a ^” R a 9 
is negligible regardless of whether 2/3-CL or full-CL para- 
meterization was employed. 

As contrasted with the Winter (16 January 1974) studied 
by Spaeth, no change in sign in R g ^ ~ R^ and/or in R^ - R^ 

occurs for <{) > 40N. The reason for the consistent sign of 

- 2.-1 

the difference R - - R 0 - .03 cal cm min which exists at 

si s 2 

virtually all latitudes in April is that the "summer like- 
effect" of an increase in effective insolation with decreased 
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Figure 14. Aadiational balance at the tropopause (R g ),at 
the ocean surface (R) , and in the tropospheric column (R a ) . 
Solid lines denote computations made with the 2/3-CL (sub- 
scripted "l") cloud model and dashed lines correspond to the 
full-CL model, (subscripted "2"). All computations are for 
16 April 1974 in ly min -1 x 10 2 . 



no 



cloud-cover is dominant at all <p £ 65N at this time of the 

year (when the declination is 8.5 degrees), whereas there 

is a smaller retention of IR net flux by the troposphere when 

the cloud-cover is reduced to 2/3-CL. These effects apply at 

all latitudes under study here, but undoubtedly, there could 

be a change of sign in both R , - R „ and R, - R_ in lati- 

sl s2 1 2 

tudes further removed from the subsolar latitude for 16 April 
1974 . 

C. CROSS-SEASONAL EFFECTS 

Tables X, XI and XII show the cross-seasonal effects in 
R , R and R centered on 16 April 1974. It is clear that 

S 3 . 

increases at each <p in the Northern Hemisphere in the time 

sequence 16 January through 16 July 1974. A similar effect 

applies to R £ (full-CL case) through the same time sequence. 

Similar results apply to the time sequence for the sizable 

increases of net flux (R) at the earth’s surface, both for 

the 2/3-CL and full-CL cases. These cross-seasonal results 

tend to confirm the validity of the radiative model. 

In the Southern Hemisphere latitudes of Tables X, XI and 

XII, there is a systematic trend toward decreasing R -values 

s 

in the time sequence 16 January through 16 July 1974. The 
same trend applies to R-values across the time sequence. 
Considering that this sequence corresponds to the periods of 
s umme r- f a 1 1-win t e r in the Southern Hemisphere, this cross- 
seasonal trend also appears to afford qualitative support to 
the radiational model tested here. 
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D. COMPARISON OF NET FLUX AT THE TOP OF THE MODEL 

ATMOSPHERE WITH SATELLITE OBSERVATIONS 

Based upon observations of Raschke et al (1973) of 
F2(RAS) and of the albedo ALB(RAS) over the latitude-range 
of this study, the net flux parameter at the top of the 
atmosphere was defined for the NIMBUS III data-period (1-15 
May 1969) by Raschke et al (1973), as 

RN = QAVE* -p"!! ( 1-ALB (RAS ) ) - F2(RAS) (7-3) 

Raschke's computed values of the RN for this data period 
used a slightly different solar constant than that employed 
in this study (S = 2.00 ly (min) ^ ) . Hence for purposes of 
consistency in comparisons, our average extraterrestrial in- 
solation, QAVE * ( 2 . 00 / 1 . 9 2 ) was used in the definition of RN 
in Eq . 7-3 ) . 

Values or RN were computed after Raschke at each five- 
degree latitude intersection on the four meridians where 
comparisons with the present model computations were possible. 
For the corresponding radiative model computations of this 
study, the following analog to Eq . 7-3 was adapted. 

RNMOD = QAVE* ( 1-ALBMOD ) - FF2 ( 7 - A ) 

FF2 represents a parameter analogous to F2(RAS), values 
of which were compared in Section III. E. on a latitudinal 
basis. Likewise, albedo values were compared on a lati- 
tudinal basis in Section IV. G. 
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The computational values of RNMOD for 16 April 1974 were 
computed using both the 2/3-CL and the full-CL cloud para- 
met erizat ions • After obtaining averages across the four 
meridians, the latitudinal comparisons are depicted in Table 
XIII . 

The first three columns of Table XIII show that the lati- 
tudinal distribution of RNMOD is somewhat smaller than the 
values of RN (after Raschke) at all latitudes except 55N. 

The difference 



RN - RNMOD 

is generally positive because of the consistently small albedos 

of the Raschke observations. This difference is considerably 

larger for the full-CL case, when the model albedo was 

recognized as being excessively high* 

A comparison of the mean radiative balance R g at the 

tropopause k=2 has been made with that at k=0, which has been 

denoted RNMOD. The results have been tabulated as a function 

of latitude in Table XIII. The R -values for the earth- 

s 

troposphere have been extracted from Table X and appear in 
columns four and five of Table XIII. Recall that 

R = QAVE (1-ALBMOD) - F * . (7-5) 

s l 

In comparing Eqs. 7-4 and 7-5, it is clear that the incoming 
insolation used in the computation of R^ has been reduced by 
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4% of an amount which is approximately 



QAVE - QREF = .3602 ly (min) -1 (7-6) 

averaged over all latitudes. However, exceeds RNMOD(2/3) 

by .017 ly (min) ^ which must therefore be accounted for by 
the fact 



FF2 - F 2 * = .023 - .006 + . 04 (QAVE-QREF) (7-7) 

for the 2/3-CL case. Approximately the same difference 
exists between comparisons for the full-CL cases. 

The IR-flux divergence of (7-7) is due to the presence 
of water-vapor and CO ^ emitter masses operating in conformity 
with the F 2 * model (Eq. 3-5). These masses have not been 
subjected to compensating solar insolation absorption in the 
stratosphere, and the 4% solar absorption by oxygen and ozone 
is not sufficient to establish a stratospheric balance. 
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Mid-Seasonal Dates 



Lat . 


16 Jan 


1974 


16 April 1974 


16 July 


1974 




R i 

s 1 


R 9 

s2 


R i 

si 


R 0 
s 2 


R -i 

si 


R 9 

s2 


2 OS 


.0681 


.0216 


-.0229 


-.0469 


-.1012 


-.1124 


15 


.0552 


.0052 


-.0119 


-.0418 


-.0817 


-.1010 


10 


.0488 


-.0039 


. 0002 


-.0358 


-.0546 


-.0772 


5 


.0475 


.0046 


. 0171 


-.0235 


-.0428 


-.0616 


0 


.0182 


-.0270 


. 0353 


-.0079 


. 0132 


-.0171 


5 


-.0050 


-.0478 


. 0283 


. 0243 


. 0297 


-.0062 


10 


-.0140 


-.0428 


. 0198 


-.0364 


. 0521 


.0145 


15 


-.0489 


-.0736 


.0225 


-.0335 


.0635 


.0211 


20 


-.0813 


-.1010 


.0204 


-.0368 


.0627 


.0136 


25 


-.0938 


- . 1058 


. 0256 


-.0126 


. 0825 


.0380 


30 


-.0998 


-.1041 


. 0547 


. 0293 


.0945 


.0560 


35 


-.1497 


-.1551 


. 0488 


.0301 


. 1018 


.0607 


AO 


-.1689 


-.1706 


.0540 


.0296 


.1254 


.0888 


45 


-.1884 


-.1835 


. 0427 


.0096 


. 0884 


.0612 


50 


-.2062 


-.1969 


-.0055 


-.0362 


. 1461 


. 1038 


55 


-.1798 


-.1671 


. 0016 


-.0146 


.1063 


.0766 


60 


-.1932 


-.1759 


-.0428 


-.0591 


. 1063 


. 1315 


65N 


- .2265 


-.2265 


-.0089 


-.0374 


. 1990 


.1640 


Wt. Avg. 


- . 1000 


-.1139 


. 0227 


-.0153 


.0523 


.0227 



Table X. Cro ss- se as ona 1 comparison of zonally-averaged 
values of R as a function of latitude. The 

g 

averaged values for the 2/3-CL case are denoted with a 
subscripted 1 and the full-CL case values by a subscripted 
2. (All values in ly (min) . 
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Mid-Seasonal Dates 



Lat . 


16 Jan 


1974 


16 April 


1974 


16 July 


1974 




R, 


R„ 


R, 


R 


R, 


R 




1 


2 


1 


2 


1 


2 


2 OS 


. 1622 


.1123 


.0813 


.0514 


.0598 


.0546 


15 


. 1624 


.1125 


. 1015 


.0681 


.0954 


.0852 


10 


. 1780 


. 1310 


.1299 


.0957 


. 1276 


. 1134 


5 


. 1902 


. 1537 


. 1438 


. 1049 


.1375 


. 1224 


0 


. 1802 


.1443 


. 1164 


.1376 


.1732 


.1486 


5 


.1697 


.1387 


. 1858 


. 1294 


. 1773 


.1463 


10 


.1645 


.1454 


. 1911 


.1460 


. 2041 


.1726 


15 


.1374 


.1226 


.1933 


.1470 


.2216 


.1866 


20 


.1123 


.1022 


.2018 


.1569 


. 2373 


.1987 


25 


.0807 


.0738 


.2079 


.1769 


. 2579 


. 2225 


30 


.0475 


.0429 


.2153 


.1915 


.2616 


.2285 


35 


.0152 


.0145 


. 1889 


.1684 


. 2774 


.2429 


AO 


-.0092 


- o 0 0 7 2 


. 1824 


.1549 


.2927 


.2605 


45 


-.0094 


-.0068 


. 1670 


.1343 


. 2364 


. 2102 


50 


-.0275 


-.0200 


.1256 


o 09 52 


. 2747 


.2309 


55 


-.0409 


-.0300 


. 1062 


. 0850 


. 1816 


.1438 


60 


-.0868 


.0683 


.0818 


.0603 


. 2725 


.2359 


65N 


-.0868 


-.0868 


. 0614 


. 0206 


. 2637 


.2154 




Wt . Avg . 


. 0668 


.0581 


. 1684 


.1333 


. 1879 


. 1349 



Table XI, Cross-seasonal comparison of zonally-averaged 
values of R as a function of latitude <f) . The 
averaged values for the 2/3-CL case are denoted with a sub- 
scripted 1 and the full-CL case values by a subscripted 2. 
(All values in ly (min)~l). 
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Mid-Seasonal Dates 



Lat . 


16 Jan 


1974 


16 April 1974 


16 July 


1974 




R i 

a 1 


R a 2 


R i 

a 1 


R a 2 


R i 

al 


R 0 
a2 


2 0 S 


-.0942 


-.0913 


-.1042 


-.0983 


-.1611 


-.1674 


15 


-.1073 


-.1075 


-.1134 


-.1099 


-.1772 


-.1863 


10 


-.1292 


-.1347 


-. 1297 


-.1315 


-.1824 


-.1907 


5 


-.1426 


-.1492 


-.1267 


-.1284 


-.1804 


-.1842 


0 


-.1621 


-.1720 


-.0811 


-. .455 


-.1603 


-.1659 


5 


-.1746 


-.1864 


-.1575 


-.1671 


-.1476 


-.1526 


10 


-.1787 


-.1882 


-.1713 


-.1824 


-.1521 


-. 1582 


15 


-.1863 


-.1962 


-.1708 


-.1805 


-.1582 


-. 1656 


20 


-.1937 


-.2031 


-.1814 


-.1937 


-.1748 


-.1852 


25 


-.1745 


-.1796 


-.1823 


-.1895 


-.1756 


-. 1847 


30 


-.1472 


-.1471 


-.1606 


-.1622 


-.1672 


-.1727 


35 


-.1648 


-.1696 


-.1401 


-. 1383 


-.1757 


-. 1824 


40 


-.1597 


-.1634 


-.1284 


-.1253 


-.1657 


-.1718 


45 


-.1791 


-.1768 


-.1243 


-.1247 


-.1482 


-.1491 


50 


-.1788 


-.1769 


-.1311 


-.1314 


-. 1286 


-.1274 


55 


-.1065 


-.1371 


-.1046 


-.0996 


-.0754 


-.0675 


60 


-.1397 


-.1077 


-.1246 


-. 1194 


-.1123 


-.1045 


65N 


-.1668 


-.1397 


-.0703 


-.0580 


-.0650 


-.0517 


Wt . Avg . 


-.1668 


-.1720 


-.1457 


-.1486 


-.1626 


-.1577 



Table XII. Cross-seasonal comparison of the zonally-averaged 
values of R as a function of latitude <p . The 
averaged values for the 2/3-CL case are denoted with a sub- 
scripted 1 and the full-CL case values by a subcripted 2. 

(All values in ly (min) - -*-). 
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La t . 


RNMOD 

2/3-CL Full-CL 


Ocean- 

R i 

s 1 


Troposphere 

R s2 


RNRAS 


20 S 


-.0523 


-.0762 


-.0229 


-.0469 


. 0369 


15 


-.0394 


-.0700 


-.0119 


-.0418 


.0554 


10 


-.0227 


-.0588 


. 0002 


-.0358 


.0745 


5 


.0076 


-.0272 


.0171 


-.0235 


. 1083 


0 


.0093 


-.0377 


.0353 


-.0079 


.1304 


5 


.0134 


-.0392 


.0283 


. 0243 


. 1545 


10 


.0115 


-.0448 


.0198 


-.0364 


.1517 


15 


.0172 


-.0390 


.0225 


-.0335 


.1337 


20 


.0161 


-.0412 


.0204 


-.0368 


.1272 


25 


.0232 


-.0150 


. 0256 


-.0126 


.1150 


30 


.0467 


.0213 


.0547 


.0293 


.1059 


35 


.0380 


.0193 


.0488 


.0301 


.0779 


40 


.0429 


.0185 


.0540 


.0296 


.0770 


45 


.0296 


-.0036 


.0427 


.0096 


.0570 


50 


-.0122 


-.0430 


-.0055 


-.0362 


.0251 


55 


-.0152 


-.0314 


. 0016 


-.0146 


-.0216 


60 


-.0340 


-.0505 


-.0428 


-.0591 


-.0489 


65N 


-.0483 


-.0770 


-.0089 


-.0374 


-.0946 


Wt . Avg . 


.0062 


-.0328 


. 0227 


-.0153 


.0818 



Table XIII, Comparison of the Ocean-Troposphere Net Radia- 
tion R at k=2, with the model computed Net 
Radiation at the top of the atmosphere RNMOD and with 
climatological Net Radiation at the top of the atmosphere 
from Raschke et al (1973) RNRAS. R and RNMOD have been 
computed using both the 2/3-CL and full-CL par amet eriza- 

tions. R - is from 2/3-CL model and R 0 is from the full-CL 
, ^ s 1 s 2 
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VIII 



ZONALLY-AVERAGED TROPOSPHERIC 
AND OCEANIC HEAT BUDGETS 



FOR 16 APRIL 1974 

A. THE TROPOSPHERIC HEAT BUDGET 

By averaging over the tropospheric columns for the four 

meridians displayed in Figs, 7, 8, 9 and 10 latitude by 

latitude, the latitudinal distribution of E + Hp results were 

obtained and are presented in Fig. 15. R , the tropospheric 

a 

radiative net cooling is listed in Table XII. The zonally- 
averaged E + Hp values computed from Eqs. 5-7 and 5-13 are 
displayed as a function of latitude in Table XIV Q The re- 
sulting heating rate of the troposphere may then be expressed 
as a function of <p by the right side of Eq. 8-1, 





Qva 


+ S = 
a 


R + 
a 


(E + Hp) 


• 






(8- 


1) 


La t . 


-20. 


-15. 


-10. 


-5. 


0.0 


5. 


10. 


15. 


20. 


E+Hp 


.2019 


. 1059 


.0815 


.0167 


. 0430 


.0355 


.0855 


. 1758 


.0645 


La t . 


25. 


30. 


35. 


40. 


45. 


50. 


55. 


60. 


65. N 


E+H r 


.0849 


.0439 


. 1499 


. 1852 


.1966 


.1318 


.1957 


.0642 


. 1314 



Weighted Average = .1032 

Table XIV. Zona lly-averaged values of E+Hp as a function 
of latitude. (All values in ly (min) -1 ) . 

is the storage heating rate of the column, and is the 

required flux-divergence of heat compatible with the heat 
balance equation for the column, for the observed values of 
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, E+Hp and S^. Figure 15 depicts each quantity (R^, 

(E + H p ) and (Q + S ) ) as a function of latitude for the 

1 va a 

two cloud-model cases. 

As noted in Section VII. B., - R £ over the range 

of latitudes in this study, with little variation in R as 

a 

a function of latitude. The distribution of E + Hp as a 
function of latitude is independent of cloud-model. There- 
fore, Fig. 15 shows only one curve for E + Hp . The curves 
of Q + S are also shown in Fig. 15, as the residual of 
the right sides of Eq. 8-1 for both the 2/3-CL and the full- 
CL cases, respectively. Note that (Q va + ^^htly 

exceeds (Q + S ) n from 10S to 30N as a result of the cor- 
va a 2 

responding excess of R^ over R^ f° r the same latitude 

range. Note that both (Q + S ) distributions reach maxi- 

va a 

mum positive values at <f> = 50N, with a general positive 
trend from 35N to 55N. This feature follows from Eq. 8-1, 
since (Q va + s a ) must follow the trend of (E + Hp ) where 
the latitudinal variation of the latter function peaks more 
sharply, relative to the more or less uniform behavior of 

R (<f>) . 

a 

The cosine-weighted value of Q + S = -.0425 ly min 

va a 

2 

per cm column. has not been computed in this study; 

however, if the Q 4* S is attributed to mean storage- 
* va a 

cooling S of the troposphere alone, the daily storage rate 
a 

corresponds to a temperature-change rate given by 

R~ + (E + H ) 

(|~) = 4.1 — x 1440 °C/ day (8-2) 
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Figure 15* Distribution of R a , E+H r and of Q va + ^ a as a 
function of latitude for the two cloud-model cases. 

(solid curves represent the 2/3-CL case and dashed curves 
represent the full-CL case) 
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with AP mb = 800mb in the troposphere. The resultant cooling 
rate (8-2) over the tropospheric depth of 800mb, considering 
that zero lateral flux-divergence applies, is approximately 

(§“) = -.31°C (day)" 1 

2 

averaged over the mean cm tropospheric column, 

B. THE LATITUDINALLY-AVERAGED HEAT BUDGET OF THE OCEAN 
The latitudinally-averaged heat-budget of the ocean 
water-mass obeys the following thermodynamic relationship 
for both the 2/3-CL and full-CL cases, 

Q + S = R + (- (E + H r ) ) (8-3) 

vo o r 

Note that Eq. 8-3 is the analog of Eq. 8-1 for the atmosphere, 
with -(E + Hp) representing the average turbulent cooling 
rate for the water-mass by mixing across the sea-air inter- 
face. 

The terms of the right side of Eq, 8-3 have been computed 
for each latitude and averaged zonally, with R and E + Hp 
shown in Table XI and XI V , r es p e c t i ve ly , 4- has been 

computed as the residual of these two terms as a function of 
latitude and the entire set of terms in Eq. 8-3 are graphed 
in Fig. 16 for both cloud cases. 

The primary ocean-mass heating-function is R, the radia- 
tive heating rate at the ocean-air interface. As discussed 
in Section VII, Rp > R^ for all latitudes considered in this 
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study. Therefore, according to Eq. 8-3, (Q yo + S q ) must 

be greater than (Q + S ) n by the same amounts since 

vo o 2 y 

-(E + = + ^ 2 * This means that the internal heat- 

ing rate of the ocean-mass, as given in Eq. 8-3, is greater 
in the 2/3-CL case than in the full-CL case for all latitudes 
considered in this study. is the required mean oceanic 

heat-flux divergence at latitude <p . Fig. 16 supports the 
conclusion of Von der Haar and Oort (1973) which requires 
increased lateral transport-divergences in the oceans, 

based upon the recently observed satellite values of the 
global albedo. The values for (Q yo + S^), from this study 
seem to corroborate the satellite findings. 

The mean value of (Q + S ) computed between 20S-65N 

vo o 1 r 

is 

(Q + S ) _ = .0652 ly min" 1 
x vo o 1 J 

which corresponds to a mean heating rate in the water-mass 
column. It should be noted that this heating rate is 
noticeably larger, and of opposite sign, than the corres- 
ponding atmospheric effect. 



12 3 



/ 







Figure 16. Distribution of R,-(E+H r ) and of Q +S as a 
function of latitude for the tv/o cloud-model cases. (solid 
curves represent the 2/3-CL case and dashed curves represent 
the full-CL case) 
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IX 



CONCLUSIONS 



The present work continues the use of the radiational 
two-layer cloud model adapted from Spaeth's mid-winter study 
to FNWC gridpoint data valid over the four meridians for 
16 April 1974. The use of the 2/3-CL parameterization, as 
compared to the full-CL after Smagorinsky (1968), for the 
specification of cloud-layer amounts gave the best verifi- 
cation for the computed global albedo as a function of lati- 
tude when compared with observations from NIMBUS III 
satellite climatology comparable to that of 16 April 1974. 
The use of the 2/3-CL parameterization also gave good veri- 
fication of computed long-wave flux to space when compared 
to the latitudinal distribution of the same element as given 
by satellite climatology. 

The chief shortcoming of the computed global albedo 
versus the satellite measured albedo occurred in tropical 
latitudes (20S-20N). Here, it appears that cloud-elements 
do not behave as reflecting sheets for insolation but rather 
as focusing devices through the intra-cloud spaces. The use 

of 50-60% model cloud reflectances seems excessive in these 

/ 

latitudes. As a practical alternative it may be reasonable 
to tune the cloud reflectances until agreement of the com- 
puted global albedos with satellite measurements is 
achieved in the latitude range 20S-20N. 
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In the topics, all analyses are subject to da ta-spar s i ty 
and hence are largely dependent upon climatology and a rela- 
tively few radiosonde reporting stations. An apparent weak- 
ness in this study is that of overspecifying cloud-amounts 
through the use of Smagorinsky type formulas. This could be 
a result of inadequately resolved relative humidities, which 
in general might well be lower than inferred by the FNWC 
analyses • 
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